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. Lý thuyết cơ bản về quy hoạch tuyến tính 
. Giải thuật đơn hình 

. Bài toán đối ngẫu 

. Ứng dụng quy hoạch tuyến tính 

. Bài tập tổng hợp 


Lý thuyết cơ bản về quy hoạch tuyến tính 

Chương này trình bày cách xây dựng mô hình quy hoạch tuyến tính của 
những bài toán dạng đơn giản. Đây là những kiến thức quan trọng để xây 
dựng mô hình cho những bài toán phức tạp hơn trong thực tẾ sau này. Các 
khái niệm về “? lồi? đuợc trình bày để làm cơ sở cho phương pháp hình 
học giải quy hoạch tuyến tính. Một ví dụ mở đầu được trình bày một 
cách trực quan để làm rõ khái niệm về phương án tối ưu của quy hoạch 
tuyến tính. 


LÝ THUYẾT CƠ BẢN VỀ QUY HOẠCH TUYẾN TÍNH 


GIỚI THIỆU BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 

Có thể tạm định nghĩa quy hoạch tuyến tính là lĩnh vực toán học nghiên 
cứu các bài toán tối ưu mà hàm mục tiêu (vấn đề được quan tâm) và các 
ràng buộc (điều kiện của bài toán) đều là hàm và các phương trình hoặc 
bất phương trình tuyến tính. Đây chỉ là một định nghĩa mơ hồ, bài toán 
quy hoạch tuyến tính sẽ được xác định rõ ràng hơn thông qua các ví dụ . 


Các bước nghiên cứu và ứng dụng một bài toán quy hoạch tuyến tính 
điển hình là như sau : 


a- Xác định vấn đề cần giải quyết, thu thập dữ liệu. 
b- Lập mô hình toán học. 


c- Xây dựng các thuật toán để giải bài toán đã mô hình hoá bằng ngôn 
ngữ thuận lợi cho việc lập trình cho máy tính. 


d- Tính toán thử và điều chỉnh mô hình nếu cần. 


e- Áp dụng giải các bài toán thực tế. 


Bài toán vốn đầu tư 


Người ta cần có một lượng (tối thiểu) chất dinh dưỡng i¡=1,2,..,m do các 
thức ăn j=1,2,...,n cung cấp. Giả sử : 


aij là số lượng chất dinh dưỡng loại ¡ có trong 1 đơn vị thức ăn loại j 
(¡=1,2,...,m) và (j=1,2,..., n) 
bi là nhu cầu tối thiểu về loại dinh dưỡng ¡ 
cj là giá mua một đơn vị thức ăn loại j 
Vấn đề đặt ra là phải mua các loại thức ăn như thế nào để tổng chỉ phí 
bỏ ra ít nhất mà vẫn đáp ứng được yêu cầu về dinh dưỡng. Vấn đề được 
giải quyết theo mô hình sau đây : 
Gọi xj 0 (j= 1,2,...,n) là số lượng thức ăn thứ j cần mua . 
Tổng chi phí cho việc mua thức ăn là : 
†+ 
2= » G77 S 6l00iT-.6a0a SE ¿uê TT nh 


Vì chi phí bỏ ra để mua thức ăn phải là thấp nhất nên yêu cầu cần được 
thỏa mãn là : 


min z — : 6723 = Cíđ 1? ri ba19 de + Cx„#„ 

Lượng dinh dưỡng ¡ thu được từ thức ăn 1 là : ai1x1 (¡=1 m) 
Lượng dinh dưỡng ¡ thu được từ thức ăn 2 là : ai2x2 

Lượng dinh dưỡng ¡ thu được từ thức ăn n là : ainxn 

Vậy lượng dinh dưỡng thứ ¡ thu được từ các loại thức ăn là : 


ai†1x1+al2x2+...+ainxn (I=1 m) 


Vì lượng dinh dưỡng thứ ¡ thu được phải thỏa yêu cầu bi về dinh dưỡng 
loại đó nên ta có ràng buỘc sau : 


ailx1+al2x2+...+ainxn bị (=1 m) 

Khi đó theo yêu cầu của bài toán ta có mô hình toán sau đây : 
° (2 

min z = 3Í G732 '0181°5EoCgiba7242gsxe + Can 


0111 + 01282 +‹:. + .01n#a > 1 
đ211 -Ƒ đạa2#2 -... + đạn#a > Ùa 
Œ(m1#1 + đm222 -E...  OmnZ#n So Dịn 
*j > 0(7 — li XS ¡Ì) 


tít 


Bài toán lập kế hoạch sản xuất 

Từ m loại nguyên liệu hiện có người ta muốn sản xuất n loại sản phẩm 
Giả sử : 

aij là lượng nguyên liệu loại ¡ dùng để sản xuất 1 sản phẩm loại j 
(¡=1,2,...,m) và (j=1,2,..., n) 

bi là số lượng nguyên liệu loại ¡ hiện có 

cj là lợi nhuận thu được từ việc bán một đơn vị sản phẩm loại j 


Vấn đề đặt ra là phải sản xuất mỗi loại sản phẩm là bao nhiêu sao cho 
tổng lợi nhuận thu được từ việc bán các sản phẩm lớn nhất trong điều 


kiện nguyên liệu hiện có. 

Gọi xj 0 là số lượng sản phẩm thứ j sẽ sản xuất (j=1,2,...,n) 

Tổng lợi nhuận thu được từ việc bán các sản phẩm là : 

CÁ 3i 78, = Ciữi:(089/đ suau + Cna2a 

Vì yêu cầu lợi nhuận thu được cao nhất nên ta cần có : 

max Z — gi Cj#j — C181 + Ca#a +...... + CnZ„ 

Lượng nguyên liệu thứ ¡=1 m dùng để sản xuất sản phẩm thứ 1 là ai1x1 
Lượng nguyên liệu thứ ¡=1 m dùng để sản xuất sản phẩm thứ 2 là ai2x2 
Lượng nguyên liệu thứ ¡=1 m dùng để sản xuất sản phẩm thứ n là ainxn 
Vậy lượng nguyên liệu thứ ¡ dùng để sản xuất là các sản phẩm là 
aI1x1+al2x2+...+ainxn 


Vì lượng nguyên liệu thứ ¡=1 m dùng để sản xuất các loại sản phẩm 
không thể vượt quá lượng được cung cấp là bi nên : 


ailx1+al2x2+...+ainxn bi (i=1,2,...,m) 
Vậy theo yêu cầu của bài toán ta có mô hình sau đây : 


†\ 
max Z = » Gi9/= 0i01cÐ E802 sua + CnZ„ 


đ11#1 + đ1a#a +...  địn#„ < bị 
0211 + 0ạa2#a --... + đan#„ ba 
(m11 + („22 Si ng (Ñmn#„ẹ X Ủm; 
%j Bo 0(7 —= D2 1x1} 


tt 


Bài toán vận tải 


Người ta cần vận chuyển hàng hoá từ m kho đến n cửa hàng bán lẻ. 
Lượng hàng hoá ở kho ¡ là sỉ (¡=1,2,...,m) và nhu cầu hàng hoá của cửa 
hàng j là dị (j=1,2,...,n). Cước vận chuyển một đơn vị hàng hoá từ kho ï 
đến của hàng j là cij 0 đồng. 


Giả sử rằng tổng hàng hoá có ở các kho và tổng nhu cầu hàng hoá ở các 
cửa hàng là bằng nhau, tức là : 


†n †+ 
3 ¡=1 8= » =í đ; 
Bài toán đặt ra là lập kế hoạch vận chuyển để tiền cước là nhỏ nhất, với 
điều kiện là mỗi cửa hàng đều nhận đủ hàng và mỗi kho đều trao hết 


hàng. 


Gọi xij 0 là lượng hàng hoá phải vận chuyển từ kho ¡ đến cửa hàng j. 
Cước vận chuyển chuyển hàng hoá ¡ đến tất cả các kho j là : 


†+ 
3271 Củ 
Cước vận chuyển tất cả hàng hoá đến tất cả kho sẽ là : 


Zz= 3-1 —1 Cụ#u 


Theo yêu cầu của bài toán ta có mô hình toán sau đây : 


min2=Š)-; 3; ¡6g 
` — Tj " đ;(j == Ì2e2 271) 
#Tj Đo 0( " 1,2,.x„t](7 = 11:28) 


QUY HOẠCH TUYẾN TÍNH TỔNG QUÁT VÀ CHÍNH TẮC 


Quy hoạch tuyến tính tổng quát 


Tổng quát những bài toán quy hoạch tuyến tính cụ thể trên, một bài toán 
quy hoạch tuyến tính là một mô hình toán tìm cực tiểu (min) hoặc cực 
đại (max) của hàm mục tiêu tuyến tính với các ràng buộc là bất đẳng 
thức và đẳng thức tuyến tính. Dạng tổng quát của một bài toán quy 
hoạch tuyến tính là : 


min/maxz = Sàn, Gi00 1) 
` G¡j# j — b;(¡ = T1) 
S3 G2; sế b;(¡ = T;)() 
337—10j#¿ > b¿(¡ € Tạ) 
#¡ > UJLJ5C ti 


z; <07€ J; (HT) 
œ;tùy ý j € J 


TìO 


tÍ 
tí 


Trong đó : 


() Hàm mục tiêu 


Là một tổ hợp tuyến tính của các biến số, biểu thị một đại lượng nào đó 
mà ta cần phải quan tâm của bài toán. 


(H) Các ràng buộc của bài toán 


Là các phương trình hoặc bất phương trình tuyến tính n biến số, sinh ra 


từ điều kiện của bài toán. 


(II) Các các hạn chế về dấu của các biến số 


Người ta cũng thường trình bày bài toán quy hoạch tuyến tính dưới dạng 
ma trận như sau : 


+1 
T2 


#n 


rìgh, 


đ11312:..31n C1 
21322...À2n C2 


r¡ìgh righ 
[II Dị 
A=lajl= bạ 
Đmm 

rìgh 


Gọi ai (¡=1 m) là dòng thỨ ¡ của ma trận A, ta có : 


mmìn/maxz(#) = cTfz(1) 
gz#.—= b;(¡ = Tị) 
LH bô b¿(¡ = 7;)(H) 
d¿œ > b¡(ñ € ]ạ) 


U07 €I 
z¡<07€ 7; (II) 
z;tùy ý j € ở 


TìO 


tí 


Người ta gỌi : 
- A là ma trận hệ số các ràng buộc. 
- c là vectơ chỉ phí (cT là chuyển vị của c) 


- b là vectƠ giới hạn các ràng buộc. 


Quy hoạch tuyến tính dạng chính tắc 


Bài toán quy hoạch tuyến tính chính tắc là bài toán quy hoạch tuyến tính 
mà trong đó các ràng buộc chỉ có dấu = và các biển số đều không âm. 
min/maxz = Xung #;#7(Ì) 


35-1 đu =b,0= 1/2,2.ai)00) 
#; > 0(7 = 1,2,...n)(TT) 


(m n) 
{ 
min/maxz(z) = cfz(1) 
Ax= 0(T) 
z > 0(1) 
rang(A)=m 


Người ta có thể biến đổi bài toán quy hoạch tuyến tính dạng tổng quát 
thành bài toán quy hoạch tuyến tính dạng chính tắc nhờ các quy tắc sau 
đầy : 


- Nếu gặp ràng buộc ¡ có dạng thì người ta cộng thêm vào vế trái của 
ràng buộc một biến phụ xn+i 0 để được dấu =. 


- Nếu gặp ràng buộc ¡ có dạng thì người ta trừ vào vế trái của ràng 
buộc một biến phụ xn+i 0 để được dấu =. 


Các biến phụ chỉ là những đại lượng giúp ta biến các ràng buộc dạng bất 
đẳng thức thành đẳng thức, nó phải không ảnh hƯởng gì đến hàm mục 
tiêu nên không xuất hiện trong hàm mục tiêu. 


- Nếu biến xj 0 thì ta đặt xj = -x'j với x'j 0 rồi thay vào bài toán. 

- Nêu biên xj là tuỳ ý thì ta đặt #; = #; — #; vỚi #;, t2, đếu 0rỗi 
thay vào bài toán. 

- Trong trường hỢp trong sỐ các ràng buộc có dòng mà vế phải của dòng 
đó là giá trị âm thì đổi dấu cả hai vế để được vế phải là một giá trị 
không âm. 

Dựa vào các phép biến đổi trên mà người ta có thể nói rằng bài toán quy 


hoạch tuyến tính chính tắc là bài toán quy hoạch tuyến tính mà trong đó 
các ràng buộc chỉ có dấu =, về phải và các biến số đều không âm. 


Ví dụ : 


Biến đổi bài toán quy hoạch tuyến tính sau đây về dạng chính tắc : 


mìnZz(#)= 2X:— ø2+:2Xs + #a4— 2X 
#1 — 23Xa + #a + 2xa + z;s < 7 
øa + 2Xa + za > —] 
2xs + #4 + 3x; > 10 
#1 + #a — 2xs + z4 = 20 


810g < +Ù 
#4 X 0 
#a,zstùy ý 


TuO 


tÍ 
TẠI 


Bằng các thay thế : 


ta đƯỢc : 


minz(#) = 2xị — (z› = #2) _ 2(z; =— z3) ` ta — 2%; 
#1 — 2(1¿ — #s ) + (8s — #x ) — 2x + #s + øe = 7 


(a — #¿ ) + 2# — #¿ ) + #¿T— #r = —] 


%( 
23 — #4) — 84 + 3g — #g = 10 


"" 2(za =1... 


2 22 là Lê 


11;45;X6;17;⁄s;12;Ÿ2 ;⁄3;X3 ;T4 Bàn 0 


tí 


hay : 


minz(#) = 2xị — (z› — #2) + 2(z; =— z3) — tụ — 2%; 
#1 — 2(Z¿ — øs ) + (8s — #x ) — 2x + #g + øạ = 7 


X8 52219092 0NG 
2(#s — #s) — 84 + 3M — 8g =10 


#1 + (%s — #s) — 2(Zs )—g) — øạ = 20 


Lộ bến) ? ?`2 


11;1ã;16;⁄7;Ÿ§;12;12 ;Ý3;T3 ;Ÿa = 0 


tí 


Phương án 


Xét bài toán quy hoạch tuyến tính chính tắc : 
min/maxz(z) = cÝz 

Ax=ö 

œ>0 


(P) 


e x=[x1 x2... xn] T là một phương án của (P) khi và chỉ khi Ax = b. 

e x=[x1 x2... xn] T là một phương án khả thi của (P) khi và chỉ khi Ax 
=bvàx 0. 

« Một phƯơng án tối ưu của (P) là một phương án khả thi của (P) mà 
giá trị của hàm mục tiêu tương ứng đạt min/max. 


ĐẶC ĐIỂM CỦA TẬP HỢP CÁC PHƯƠNG ÁN 
Khái niệm lồi và các tính chất 
a- Tổ hợp lồi 


- Cho m điểm xi trong không gian Rn. Điểm x được gọi là tổ hợp lồi của 
các điểm xi nếu: “ — 2L t2” = G18) + Ga) + e + đu” 
Œ1,Gœs,....,œ„ > ơi + œa +.... +ơ„y — Ì 
- Khi x là tổ hợp lồi của hai điểm x1, x2 người ta thường viết : 

x= xi+†(1- )x2(0 1) 


Nếu 0< <1 thì x được gọi là tổ hợp lồi thật sự. 


- Đoạn thẳng 


Tập hợp tất cả các tổ tổ hợp lồi của 2 điểm bất kỳ A,B_ Rn được gọi là 
đoạn thẳng nối A và B. Ký hiệu : 


AB={x= A+(-)Bvới [01]} 
Định lý 
Tổ hợp lồ có tính chất bắc cầu. 
b- Tập hợp lồi 


Tập con S của Rn được gọi là tập hợp lồi khi S chứa toàn bộ đoạn thẳng 
nối hai điểmbất kỳ của S. 


x†-)y Š xy, I01] 
Tập hợp rỗng và tập hợp chỉ có một phần tử được xem là tập hợp lồi. 
Định lý 


Giao của một số bất kỳ các tập hợp lồi là một tập hợp lồi. 


Định lý 


Nếu S là một tập hợp lồi thì S chứa mọi tổ hợp lồi của một họ điểm bất 
kỳ trong S. 


c- Điểm cực biên của một tập hợp lồi 
Điểm x trong tập lồi S Rn được gọi là điểm cực biên nếu không thể 
biểu diễn được x dưới dạng tổ hợp lồi thật sự của hai điểm phân biệt 


của S. 


X 


4- Đa diện lồi và tập lồi đa diện 


Đa diện lồi 


Tập hợp S tất cả các tổ hợp của các điểm x1, x2,.....xm cho trước được 
gọi là đa diện lồi sinh ra bởi các điểm đó. 


Đa diện lồi là một tập hợp lồi. 

Trong đa diện lồi người ta có thể loại bỏ dần các điểm là tổ hợp của các 
điểm còn lại. Khi đó người ta thu được một hệ các điểm, giả sử là y1, 
y2,...,yp (p m) . Các điểm này chính là các điểm cực biên của đa diện lỗi, 
chúng sinh ra đa diện lồi đó. 

Số điểm cực biên của đa diện lồi là hỮu hạn. 

Siêu phẳng - Nửa không gian 

A=[aij]m.n là ma trận cấp m.n 

Ai (¡=1,2,...m) là hàng thứ ¡ của A 

Siêu phẳng trong Rn là tập các điểm x=[x1,x2,.....,xn]T thỏa 

AIx =bi 

Nửa không gian trong Rn là tập các điểm x=[x1,x2,.....,xn]T thỏa 

Alx bi 


Siêu phẳng và nửa không gian đều là các tập hợp lồi. 


Tập lồi đa diện 


Giao của một số hữu hạn các nửa không gian trong Rn được gọi là tập 
lồi đa diện. 


Tập lồi đa diện là một tập hợp lồi. 


Nếu tập lồi đa diện không rỗng và giới nội thì đó là một đa diện lồi 


Đặc điểm của tập hợp các phƯơng án 
Định lý 


Tập hợp các phương án của một quy hoạch tuyến tính là một tập lồi đa 
diện. 
Nếu tập hợp lồi đa diện này không rỖng và giới nội thì đó là một đa diện 
lồi, số điểm cực biên của nó là hỮu hạn. 
Định lý 
Tập hợp các phương án tối ưu của một quy hoạch tuyến tính là một tập 
lồi. 
Xét quy hoạch tuyến tính chính tắc 
min/maxz(z) = cfz(1) 

Ax= 0(T) 

z > 0(IT 


Giả sử A=[aij]m.n có cấp m.n,m_ n, rang(A)=m. 


Gọi Aj (j=1,2,...,n) cột thứ j của ma trận A, quy hoạch tuyến tính chính 
tắc trên có thể viết : 


min/maxz(#) = c1Z1 + ca#a +... Ƒ Cn#„ 
#z1A! + z2 A2 +... -+m„A" =b 
#:.U 


Gọi S={x=[x1,x2,..,xn]T 0/x1A1+ x2A2+...+ xnAn=b} là tập các 
phƯơng án của bài toán. 


T Ä ^ z z 
œ9 = |z1,zÿ,...,=Đ]|'_ S là một phương án khác 0. 


Định lý 


Điều kiện cần và đủ để x0 là phương án cực biên ( điểm cực biên của S) 
là các cột Aj Ứng với z?>0 là độc lập tuyến tính. 


Hệ quả 


Số phương án cực biên của một quy hoạch tuyến tính chính tắc là hữu 
hạn. Số thành phần > 0 của một phƯƠng án cực biên tối đa là bắng m. 


Khi số thành phần > 0 của một phương án cực biên bằng đúng m thì 
phƯơng án đó được gọi là một phương án cơ sỞ. 


Định lý 


Nếu tập các phương án của một quy hoạch tuyến tính chính tắc không 
rồng thì quy hoạch tuyến tính đó có ít nhất một phương án cực biên. 


z là một phương án tối ưu của quy hoạch tuyến tính. 

x1, x2 là các phương án của quy hoạch tuyến tính. 

z là tổ hợp lồi thực sự của x1, x2 

thì x1, x2 cũng là phương án tối ưu của quy hoạch tuyến tính. 
Định lý 


Nếu quy hoạch tuyến tính chính tắc có phương án tối ưu thì thì sẽ có ít 
nhất một phƯơng án cực biên là phương án tối Ưu. 


Ví dụ : xét quy hoạch tuyến tính chính tắc 
maxz(%) = 2xị + 3xạ 


4x + 2xa + #a = 5 
#1 +äxa = Ì 


T1;T2;,T3 Đài 0 


{ 


Với hệ A1 A2 ta tính được z! = [J#J —+ 0|” 


Với hệ A1 A3 ta tính được z2=[1 0 1Ị” 


— : ÉE 
Với hệ A2 A3 ta tính được z3 = I0 = “| 


Vì các thành phần của phƯƠơng án cực biên là > 0 nên ta chỉ xét x2 và x3. 
Khi đó : 


z(x2)=2.1+3.0=2 

z(x3)=2.0+3.1/3=1 

Vậy z2 =[1 0 1] là một phương án tối ưu. 
Định lý 


Điều kiện cần và đủ để một quy hoạch tuyến tính có phƯƠơng án tối ưu 
là tập các phương án không rỗng và hàm mục tiêu bị chặn. 


Định lý 

Nếu tập các phương án của một quy hoạch tuyến tính không rỗng và là 
một đa diện lồi thì quy hoạch tuyến tính đó sẽ có ít nhất một phương án 
cực biên là phương án tối Ưu. 

Phương pháp hình học 


Từ những kết quả trên người ta có cách giải một quy hoạch tuyến tính 
hai biến bằng phương pháp hình học thông qua ví dỤ sau : 


Ví dụ : xét quy hoạch tuyến tính 


maxz(z) = 3xị + 2xa 
#Ø1—#aạ> 4 
ø1 + 2xạ < 14 
5xị + 2x¿ < 30 


1;Z2 Đà 0 


tÍ 


x2x1ABCDO 

A,B,C,D,O là các điểm cực biên. Giá trị hàm mục tiêu tại đó là : 
z(A)=3.6+2.0=18 

z(B)=3.4+2.5=22 

z(C)=3.2+2.6=18 

z(D)=3.0+2.8=8 

z(O)=3.0+2.0=0 


Phương án tối ưu của bài toán đạt được tại B : x1=4 và x2=5 


MỘT VÍ DỤ MỞ ĐẦU 


Xét bài toán quy hoạch tuyến tính : 


min z(#) — -ðxị — 4xa — 3X 
2xI+äxa+za S5 
4x + z¿ + 2xa < 11 
3x + 4xa + 2xs < 8 


T1;T2;,T3 Đa 0 


tí 


Đưa bài toán về dạng chính tắc bằng cách đưa vào các biến phụ w1, w2, 
w3_ 0( làm cho các ràng buộc bất đẳng thức thành đẳng thức ) . Ta được 


min z(#) — -ðxị — 4xạ — 3X 
2xI +3xa †+ #a T11 =5 
4x1 + za + 2xa + uạ = T1] 
XI + 4xa + 2xs + ạ = 8 


11;2;3;101;102;103 Đài 0 


tí 


Thực hiện việc chuyển vế ta được bài toán ban đầu nhƯ sau : 


min z() — -ðxị — 4xa — 3X 
111 —5— 2X] — ỔXa — #3 
10a — ll — 4x] — #ạ — 2%X3 
1a —= 8— 3x1 — 4xa — 2X 


(0) 


11;2;3;101;102;103 Đụ 0 


tí 


Một phương án khả thi xuất phát ( chưa là phương án tối ưu ) của bài 
toán là : 


xi=x2=x3=0 
w1=5w2=11w3=8 
Giá trị tương ứng của hàm mục tiêu là z(x) = 0 


Người ta sẽ cải tiến phương án xuất phát này để được một phương án 
mới tốt hơn, nó làm cho giá trị của hàm mục tiêu giểm xuống. Người ta 
làm nhƯ sau : 


Vì hệ số của x1 trong hàm mục tiêu là âm và có giá trị tuyệt đối lớn nhất 
nên nếu tăng x1 từ bằng 0 lên một giá trị dương ( càng lớn càng tỐt ) và 
đồng thời vẫn giữ x2 và x3 bằng 0 thì giá trị của hàm của hàm mục tiêu 
sẽ giảm xuống. Khi đó các biến ở vế trái của bài toán (1) sẽ bị thay đổi 
theo nhưng phải thoả 0. Sự thay đổi của chúng không ảnh hưởng đến 
sự thay đổi của hàm mục tiêu. Thực hiện ý tưởng trên ta đƯỢc : 


4+ —=5—2xị¡ >0 
0ạ¿ = l11— 4xi>0 
4a —=8—3xị>Ũ 


8ù = 5. =0 


tí 


n 
IA 


lì 
— 
LA 
LA co|oo | k2|en 


t2|ct 


_ 
IA 


SuyTâ: —. x, (dòng 1 được chọn) 


tí 


Người ta chọn z¡ = š nên nhận được một phương án tốt hơn được xác 
định như sau : 


đfgsc=:XfnuSSzTÌNI IE LỆ 


Zi = ŠWa = W3 = Ÿ 

.z h ⁄ 2 ` .ÃÀA ` =. 25 
Giá trị tương Ứng của hàm mục tiêu là z(+) = —S° 
Bước tiếp theo là biến đổi bài toán (1) thành một bài toán tương đương 
bằng cách từ dòng 1 ( dòng được chọn ) tính x1 theo các biến còn lại và 
thể giá trị nhận được vào các dòng còn lại, ta được : 


: 2 l 
min z(#) = -'Š' + Š1 + 282 T— J3 
c—y. j0 1 R) L 
T1I— 3— sUỊ_— 53#2— 373 

2 


Se Q vs) Ô: `. 


qU 


11;2;Z3;1U01;,1U2;,1U03 Su 0 


tí 


Thực hiện tương tự như trên, người ta tăng x3 từ bằng 0 lên một giá trị 
dương cho phép và đồng thời vẫn giữ x2 và w1 bằng 0 thì giá trị của hàm 
của hàm mục tiêu sẽ giảm xuống. Khi đó các biến ở vế trái của bài toán 
(II) sẽ bị thay đổi theo nhưng phải thoả 0. Ta được : 


z3 <5 ( dòng 3 được chọn ) 


Khi đó người ta chọn x3=1 nên thu được một phương án tốt hơn được 
xác định nhƯ sau : 


#2 — 11 —= 10a = 


So .IWbS=vl 
Giá trị tương ứng của hàm mục tiêu là z(x)=-13 


Bước tiếp theo là biến đổi bài toán (II) thành một bài toán tương đương 
bằng cách từ dòng 3 ( dòng được chọn ) tính x3 theo các biến còn lại và 
thế giá trị nhận được vào các dòng còn lại, ta được : 
minz(z) = -l3 + 1 + 3x¿ + U3 
#1 — 2-2W1 — 2Xa + 103 
40a —= Ì+ 2Ww1 + ð%x›a 
Hứ) =l+äw + #23 —-2Wa 


(II) 


11;2;3;1U01;102;103 Mu, 0 


tí 


Đến đây vì không có hệ số nào của hàm mục tiêu là âm nên không thể 
làm giảm giá trị của hàm mục tiêu theo cách như trên nữa. Phương án thu 
được ở bước sau cùng chính là phương án tối ưu của bài toán. 

Đối với bài toán max, thay cho việc làm tăng biến có hệ số âm trong hàm 
mục tiêu người ta làm tăng biến có hệ số dương cho đến khi các hệ số 
trong hàm mục tiêu hoàn toàn âm. 


DẤU HIỆU TỐI ƯU 


Ma trận cơ sở - Phương án cƠ sỞ - Suy biến 


Xét bài toán quy hoạch tuyến tính chính tắc 


min /maxz(z) = cÝz 


J= Ở 
œ>0 


(P) 


a- Ma trận cƠ sở 

Người ta gọi cơ sở của bài toán quy hoạch tuyến tính chính tắc (P) là 
mọi ma trận B không suy biến (có ma trận nghịch đảo) mxm trích ra tỪ m 
cột của ma trận ràng buộc A. Các cột còn lại được gọi là ma trận ngoài 
cơ sở, ký hiệu là N. 

b- Phương án cơ sở - Phương án cơ sở khả thi 

B là một cơ sở của bài toán (P). 


Khi đó, bằng cách hoán vị các cột của A người ta có thể luôn luôn đặt A 
dưới dạng : 


A =[BN| 

Do đó, người ta cũng phân hoạch x và c như sau : 
xT =[xBxN | 

cT={[ cB cN ] 


Một phương án x của bài toán (P) thoả : 


Ax=b<e |BN| 
%1B 
t{N 


Phương án cƠ sở 


Người ta gọi mỘt phƯƠơng án cơ sở tương ứng với cơ sở B là mỘt 
phương án đặc biệt, nhận được băng cách cho : 


xN=0 


Khi đó xB được xác định một cách duy nhất bằng cách giải hệ phương 
trình tuyến tính bằng phương pháp Cramer : 


BxB=b xB=B-1b 

Phương án cơ sở khả thi 

Một phương án cơ sở là phương án cơ sở khả thi nếu : 
xB=B-lb 0 

CƠ sở tương ứng với mỘt phương án khả thi được gọi là cơ sở khả thi . 


Ví dụ : xét bài toán quy hoạch tuyến tính dạng chính tắc : 


min/maxz(#) = #1 — #a + #3 — #4 + #s + #6 
2xị + 2xa + #z; = 20 
—JXỊ + 4xa — 4x4 + #s — 10 
#1 + 2x + #s + äxa = 28 


z; > 0(7 = 1,2,...,6) 


tí 


Ma trận ràng buộc là 


+%1X2XäX4X5X6 
35U:U 2.10 
-340-401 
121900 


rìgh 


| 
J= 


F=nT 


Có thể chọn ba cột bất kỳ và kiểm chứng xem đó có thể là cơ sở không. 


Một cơ sở được chọn và sắp xếp lại là 


#5X6X3X4aX1X2 
100220 
010-4-34 
001512 


rìgh 


Các cột x5 x6 x3 tạo thành một ma trận cơ sở . Các biến tương ứng được 
gọi là các biến (trong) cƠ sở. 


Các cột x1 x2 x4 tạo thành một ma trận ngoài cơ sở. Các biến tương ứng 
được gọi là các biẾn ngoài cơ sở. 


Một phương án cơ sở khả thi của bài toán là : 


x1 X” xã x4 x5 x6 


c- Suy biến 


Một phương án cơ sở khả thi được gọi là suy biến nếu xB =B-1b 0 có 
những thành phần bằng 0. Sự suy biến là một hiện tượng thường xảy ra 
trong một số bài toán như bài toán vận tải, dòng dữ liệu, đường đi ngắn 
nhất....... Đây là hiện tượng khá phức tạp (có nhiều cách giải quyết sẽ 
được xét sau). Vì vậy trong những phần tiếp theo ta giả sử rằng phương 
án cƠ sở khả thi là không suy biến, tức là xB = B-1b > 0 ( dương thực sự 
). 


Dấu hiệu tối ưu 


Theo trên, khi một bài toán quy hoạch tuyến tính có phƯơng án tối Ưu thì 
tồn tại một cơ sở khả thi (tối ưu) B*, tức là phương án cơ sở x* tương 
ứng với B* là phương án tối Ưu. 

Vấn đề bây giờ là xác định một thủ tục để tìm B*. Chúng ta sẽ thấy rằng 
thủ tục đó được suy ra một cách trực tiếp tỪ việc chứng minh dấu hiệu 
tối ưu sau đây. 


Định lý 4 (dấu hiệu tối ưu) 
Xét bài toán quy hoạch tuyến tính chính tắc 
min/maxz(z) = cŸz 
Ax=ö 
œ>0 


Điều kiện cần và đủ để một phương án cơ sở khả thi x có dạng : 


zœsp—= Blb>0 
8y == Ú 


của bài toán là phương án tối Ưu là : 

ếy = cặ — cs>BN < 0 đối với bài toán max 
Zq = cặ — caB~1N > 0đỗi với bài toán min 
VỚI : 

A=[B|N] 

cT=[ cB |cN ] 

Người ta thƯờng gọi : 

cN là chi phí ngoài cƠ sở 

cB là chi phí cơ sở 

Z1 là chỉ phí trượt giảm 

c>B 1N là lượng gia giảm chỉ phí 


Chứng minh (cho bài toán max) 


Điều kiện đủ 

Giả sử x* là một phương án cơ sở khả thi với ma trận cơ sở B và thoả 
Si ¬ VẾPh, c2 TP] 

CN CƠN —CcpB `N <0 


thì cần chứng minh x* là phương án tối ưu, nghĩa là chứng minh rằng với 
mọi phương án bất kỳ của bài toán ta luôn có : 


Z0X) 2(X*) 


Xét một phương án khả thi x bất kỳ, x thoả : 
Ax=öb _ 
œ>0 


B là ma trận cơ sở của phương án cơ sở khả thi x* 


@b5 U-#zw <0 


B có ma trận nghịch đảo là B-1 


Bxp+Nxy =D 
#p>0xw>0 


B}Bxp+ B!Nxụ = B}b(E}B = T) 
#p>0xux >0 


+#b + BNxy =B!b 
#p>0xw>0 


Tính giá trị hàm mục tiêu đối với phương án x ta đƯỢc : 


Zz(x)= cÏx 


H3?) 
= T T .... T 
= léP C| R | — Cb#bB - CNTN 


=c? 1 !b— B 'Nxụy + C#N 

=cbB 1b— cT-B!NÑxN + cụ#N 

=cpB !1b+ (c-csB 1N) (Ù 

Vì x* là phương án cơ sở khả thi tương ứng với ma trận cơ sở B nên 


“TH SN 9 Sn) 


Tính giá trị hàm mục tiêu đối vơi phương án cơ bản x* ta đƯợc : 


z(x*)= clIx* 


% 
=Íe ‹#|| 5Ì =j=p+cRex 
ŸN 
Si s00 DVLdpeeU)(2) 
Từ (1) và (2) ta có : 


z(x) z(x*)vìcw —cbSB~1NÑ <0 
Vậy x* là phương án tối ưu. 
Điều kiện cần 

Œ % 


ga =0 >0 


rìgh ¬-- 3y ả ^ 
9 là phương án tÔi Ưu với ma trận cơ sở B, cần 


chứng minh rằng : ế.. = cụ — c2 1N <0. 
( ẽw là vectơ có n-m thành phần) 
Ta sẽ chứng minh điều này bằng phản chứng. 


Giả sử rằng tồn tại một thành phần cs của £;y mà cs > 0. DỰa vào cs 
người ta xây dựng một vectƠ x như sau : 


%9: 052.1 NXN 
#ỞỨ#N— ĐI, S 0 


Trong đó >0 và Is là một vectơ có (n-m) thành phần bằng 0, trừ thành 
phần thứ s bằng 1. Vậy 


Ở#N— ĐT, ., 0 
#pB =#p— B-!N0I, = B!b— B-!N0I, 
rìgh 
() 


Do B-1b_ 0 nên người ta có thể chọn >0 đủ nhỏ để xB > 0 
Vậy x được chọn như trên sẽ thoả : 
x 0(3) 


Ta kiểm chứng x thỏa ràng buộc của bài toán bằng cách tính : 


+B 


Ax=(B NỊ|“ | = Bxy + Nxy 


#{N 
=Bzp-— B!N0I, + NI, 

=BB 'b— B!NGI, + NÓI, 

=BB 'b— BB ' NóI, + NI, 

=b— NồI, + NI, 

= b(4) 

Từ (3) và (4) cho thấy x là một phương án khả thi của bài toán 


Bây giờ ta chỉ ra mâu thuẩn bằng so sánh giá trị hàm mục tiêu tại x và x* . 
Ta có : 


Zz(x)= cÏx 


+#b 
" T đ° =..‹ T 


=7) #i T 
=c{p#p— ` NxNn +CWy#N 
—. Í) T'a—1 SN 
=Cb#p —cpB *NXxN + CWy#N 


cư T T' m—1 ấh nà 
= Cb#p + Cy#yy — cgÐ *NÑxN + c{#N(vÌcwy#y —= 0) 


% 
= |cbex] lãi + cy —cpbBT}N xụ 
N 


=cTz+ cạy -cpbB"'N0I, 
=c†z + £0Ï,= c?z + 1,0 


= Z(x*) + 6;Ø > z(x*) ( vì e0 > 0) 


Vậy x* không phải là phương án tối ưu nên mâu thuẩn với giả thiết . 
Chú ý 

Qua việc chứng minh định lý dấu hiệu tối ưu ta thấy rằng từ một 
phương án cơ sở khả thi chưa tối ưu có thể tìm được các phương án khả 
thi càng lúc càng tốt hơn nhờ lặp lại nhiều lần công thức (*). Vấn đề 
được đặt là đại lượng được chọn như thế nào để nhanh chóng nhận 
được phương án tối Ưu. 

Bổ đề 


Xét bài toán quy hoạch tuyến tính chính tắc 


với B là một cơ sở khả thi nào đó và x0 là phương án cơ sở tương ứng, 
tức là 


XI Z0} =0 D10 


Xét Củ — En = c5 N, 


Nếu tồn tại một biến ngoài cơ sở xs sao cho ẽ;>0 với ẽ;là thành phần 
thứ s của e„ thì : 

a- Hoặc là người ta có thể làm tăng một cách vô hạn giá trị cỦa xs mà 
không đi ra khỏi tập hợp các phương án khả thi, và trong trường hợp này 
phương án tối ưu của bài toán không giới nội. 


b- Hoặc là người ta có thể xác định một cơ sở khả thi khác là Ö có 
phương án cơ sở khả thi z tương ứng với nó là tốt hơn, tức là : 


Z(x0) < ¿(#) 
Chứng minh 


Trong quá trình chứng minh định lý dấu hiệu tối ưu ta có phƯơng án mới 
được xác định nhƯ sau : 


%#@ =:Bl..>:0 


#B —=#p— E.NDTI,.=bE7b=b Ng 


rìgh 


Ký hiệu : 
NE—EBÌN 
N, là cột s của ý 


Di H nh 


Như vậy ta có : 


Hai trường hợp có thể xảy ra như sau : 


S 


a- Trường hợp W, < 0 
Trong trường hợp này xs có thể nhận một giá trị lớn tuỳ mà vẫn đảm 


bảo xB_ 0, nghĩa là x luôn luôn thoả 0. Khi đó như đã biết giá trị hàm 
mục tiêu tương ứng là 


lcbem] 


+B 
{N 


=cpB'b— B'NGI, + cN01, 

=c_B †1b— c_QB!NồI, + cn01, 

=z(z°) + cặ —cbB-1N 0I, 

= z(z9) + Z6], 

= z(x0) + £;8 

với 6,8 có thể lớn vô hạn thì giá trị của hàm mục tiêu là không giới nội. 
b- Trường hợp tổn tại ¡=1 m sao cho Ñ¡¿ > 0 

( N¡¿ > 0 là thành phần thứ ¡ của W,) 


Trong trường hợp này giá trị của >0 mà xs có thể nhận không thể tăng 
vô hạn vì phải đảm bảo xB>0. Giá trị lớn nhất 0 của mà xs có thể nhận 


được xác định nhƯ sau : 


p= min c—,Ñy ho 0} = 
Ns 


N 
(V¿ = 1 — mì) 


Phương án cơ sở khả thi mới có các thành phần như sau : 
+#b — b == 0N, 


ỞN— 61, 


rìgh 


và giá trị hàm mục tiêu tương ứng là : 
z() = 2z(z) + 6, > z(z°) 
Ghi chú : 


Trong trường hợp bài toán không suy biến, nếu Ø được xác định một cách 
duy nhất thì phương án mới z có đúng m thành phần khác 0. Thật vậy : 


- Biến xs đang bằng 0 trong phương án x0 trở thành dương thật sự vì 
HH r) 


- Biến xr đang dương thật sự bây giỜ nhận giá trị : 


#„ = b„ — ÔN = b„ — b: Nh = b„ — b„ = Ö 


Vậy phương án mới z là một phương án cơ sở. Nó tương Ứng với cƠ sở 
ở B được suy ra từ B bằng cách thay thế cột r bằng cột s. 


Người ta nói rằng hai cơ sở B và là kề nhau, chung tương ứng với 
những điểm cực biên kề nhau trong tập hợp lồi S các phương án khả thi 
của bài toán. 


CÂU HỎI 

1- Trình bày các bước nghiên cứu một quy hoạch tuyến tính. 

2- Định nghĩa quy hoạch tuyến tính chính tắc. 

3- Trình bày khái niệm về phương án của một quy hoạch tuyến tính. 


4- Trình bày cơ sở lý thuyết của phương pháp hình học giải một quy 
hoạch tuyến tính hai biến. 


BÀI TẬP CHƯƠNG 1 


1- Một nhà máy cán thép có thể sản xuất hai loại sản phẩm : thép tấm và 
thép cuộn. Nếu chỉ sản xuất một loại sản phẩm thì nhà máy chỉ có thể 
sản xuất 200 tấn thép tấm hoặc 140 tấn thép cuỘn trong một giỜ . Lợi 
nhuận thu được khi bán một tấn thép tấm là 25USD, một tấn thép cuộn là 
30USD. Nhà máy làm việc 40 giỜ trong một tuần và thị trường tiêu thụ 
tối đa là 6000 tấn thép tấm và 4000 tấn thép cuộn . 


Vấn đề đặt ra là nhà máy cần sản xuất mỗi loại sản phẩm là bao nhiêu 
trong một tuần để đạt lợi nhuận cao nhất. Hãy trình bày bài toán quy 
hoạch tuyển tính cho vẫn để trên. 


2- Có 3 người cùng phải đi một quảng đường dài 10km mà chỉ có một 
chiếc xe đạp một chổ ngồi. Tốc độ đi bộ của người thứ nhất là 4km/h, 
người thứ hai là 2km/h, người thứ ba là 2km/h. Tốc đỘ đi xe đạp của 
người thứ nhất là 16km/h, người thứ hai là 12km/h, người thứ ba là 
12kmh. 


Vấn đề đặt ra là làm sao để thời gian người cuối cùng đến đích là ngắn 
nhât. Hãy trình bày bài toán quy hoạch tuyến tính cho vân để trên. 


3- Một nhà máy sản xuất ba loại thịt : bò, lợn và cừu với lượng sản xuất 
mỗi ngày là 480 tấn thịt bò, 400 tấn thịt lợn, 230 tấn thịt cừu. Mỗi loại 
đều có thể bán được ở dạng tươi hoặc nấu chín. Tổng lượng các loại 
thịt có thể nấu chín để bán là 420 tấn trong giờ và 250 tấn ngoài giỜ. Lợi 
nhuận thu được từ việc bán một tấn mỗi loại thịt được cho trong bảng 
sau đây : 


Tươi Nấu chín trong giờ Nấu chín ngoài giờ 
Bò lô 14 11 
Lợn 4 12 7 
Cừu 4 13 9 


Hãy trình bày bài toán quy hoạch tuyến tính để nhà máy sản xuất đạt lợi 
nhuận cao nhất. 


4- Một xưởng mộc làm bàn và ghế. Một công nhân làm xong một cái bàn 
phải mất 2 giờ, một cái ghế phải mất 30 phút. Khách hàng thường mua 
nhiều nhất là 4 ghế kèm theo 1 bàn do đó tỷ lệ sản xuất giữa ghế và bàn 
nhiều nhất là 4:1. Giá bán một cái bàn là 135USD, một cái ghế là 50USD. 
Hãy trình bày bài toán quy hoạch tuyến tính để xưởng mộc sản xuất đạt 
doanh thu cao nhất, biết rằng xưởng có 4 công nhân đều làm việc 8 giỜ 
mỗi ngày. 


5- Một nhà máy sản xuất hai kiểu mũ. Thời gian để làm ra mỘt cái mũ 
kiểu thứ nhất nhiều gấp 2 lần thời gian làm ra một cái kiểu thứ hai. Nếu 


sản xuất toàn kiểu mũ thứ hai thì nhà máy làm được 500 cái mỗi ngày. 
Hàng ngày, thị trường tiêu thụ nhiều nhất là 150 cái mũ kiểu thứ nhất và 
200 cái kiểu thứ hai. Tiền lãi khi bán một cái mũ kiểu thứ nhất là 8USD, 
một cái mũ thứ hai là 5USD. Hãy trình bày bài toán quy hoạch tuyến tính 
để nhà máy sản xuất đạt lợi nhuận cao nhất. 


6- Trong hai tuần một con gà mái đẻ được 12 trứng hoặc ấp được 4 
trứng nở ra gà con. Sau 8 tuần thì bán tất cả gà con và trứng với giá 
0,6USD một gà và 0,1USD mội trứng. Hãy trình bày bài toán quy hoạch 
tuyến tính bố trí 100 gà mái để trứng hoặc ấp trứng sao cho doanh thu là 
nhiều nhất. 


7- Giải những bài toán quy hoạch tuyến tính sau đây bằng phương pháp 
hình học : 


1"aX2 — #1 — #2 


ImMInW — —Z1 + #2 

XI +2 >3 

—#1 — 2X <Ổ 
zø++2xa> 4 

ø1— 2xạ <4 
t1 02: e1 

—#1 +#a< Ì1 

¡<5 
07705: SS40 


8)- #3 Sð  b} 


{ tí 


maxz = 5xị †xaạ Mminw = -2XỊ — #2 


#1 — 2xạ > 2 ø1 + 2xạ <6 
=2xi + 0xa c2 Zi—#a>Š 
#1,#atuy ý 195 .2.0 
CÌề d)- 


tí lNt 


maxz — 3x — 4xa 


maxz = xị + 2X› #1—#ạ > —-4 
2xi +za <2 2xịi +#a < 14 
3xi +4xa > 1 øạ SỐ 
41,2 . 0 #1 kế 6 
e)- Ð- Z1, >0 


tí 
tít 


min/maxz(z) = 4x + 3xa 
2x1 — äxa > —12 
2xị + äxa < 24 
3X — #2 < 14 
z1 +4xa>9 


2xi+zạ54 
sâm 


®ịzZy =0 


ti 


Giải thuật đơn hình 

Chương này trình bày một cách chỉ tiết nội dung của giải thuật đơn hình. Sau phần cơ sở 
lý thuyết của giải thuật là các ví dụ tương ứng. Các ví dụ được trình bày đúng theo các 
bước của giải thuật. Kiến thức trong chương này cần thiết cho việc lập trình giải quy 
hoạch tuyến tính trên máy tính. 


GIẢI THUẬT ĐƠN HÌNH CƠ BẢN 


Chương này trình bày một phương pháp để giải bài toán quy hoạch tuyến tính đó là 
phương pháp đơn hình. Phương pháp đơn hình được George Bernard Danữig đưa ra năm 
1947 cùng lúc với việc ông khai sinh ra quy hoạch tuyến tính. Đây là một phương pháp 
thực sự có hiệu quả để giải những bài toán quy hoạch tuyến tính cở lớn trong thực tế. 
Với cách nhìn hiện đại ý tưởng của phương pháp đơn hình rất đơn giản. Có nhiều cách 
tiếp cận phương pháp đơn hình, chương này trình bày một trong các cách đó. 


CƠ sở xây dựng giải thuật đơn hình cơ bản 
Xét bài toán quy hoạch tuyến tính chính tắc : 
maxz(#) = ca 
Äx= 
>0 


Giả sử rằng B0 là một cơ sở khả thi xuất phát của bài toán ( không nhất thiết là m cột 
đầu tiên của ma trận A ) . Thuật toán đơn hình cơ bản được xây dựng dựa trên các bước 
sau : 

a-Gán B = B0 và l=0 ( số lần lặp ) 

b-l = l+1 


c-Với cơ sở hiện thời B tính : 


: phương án cơ sở khả thi tương Ứng 


b=B'b 
ấy =cwy — cyB1N: dấu hiệu tối ưu 


d-Nếu ế% = cặ — csB 1N < 0 thì giải thuật dừng và bài toán có phương án tối ưu là x 


Ngược lại, nếu tồn tại s sao cho ẽ¿ > 0 ( ẽ; là thành phần thứ s của ẽ„y) thì sang bước e 
e-Tính: A4; = B~1!A,( As là cột thứ s của A ) 
Nếu 4; < 0 thì giải thuật dừng và phương án tối ưu không giới nỘi. 


Ngược lại, nếu tỒn tại ã¡¿ € 4; mà ã¡¿ > 0thì tính : 


>0 2Ÿ si) 


z, = min gay > 0} = ý 
ø;s là các thành phần của 4;. 

z„ là thành phần thứ s của phương án mới z. 

f-Gọi xt là biến tương ứng với cột thứỨ r của cơ sở B. Khi đó biến xs sẽ nhận giá trị 

ø„ > 0 ( vào cơ sở ), biến xt sẽ nhận giá trị z¿ = 0( ra khỏi cơ sở ). Như vậy phương án 
mới z tương ứng với cơ sở mới  ( thay đổi cơ sở ) được xác định nhƯ sau : 

B=B {t}-{s} 

g-Gán B = Ö và quay về b. 


Về mặt hình học, giải thuật này được hiểu như là một quá trình duyệt qua các điểm cực 
biên của đa diện lồi S các phương án khả thi của bài toán. 


Về mặt đại số, giải thuật này được hiểu như là một quá trình xác định một chuỗi các ma 
trận cơ sở kề B0 B1 B2......... mà các phƯƠng án cƠ sở tương ứng x0 x1 x2........ là ngày 
càng tốt hơn, tức là : 


Z7 XU) <7 X1) < (X2 Ìvovb 
Chú ý : 


Nếu cơ sở ban đầu B0 chính là m cột đầu tiên cỦa ma trận A thì trong giải thuật trên t 
chính là r. 


Định lý về sự hội tụ 

Với giả thiết bài toán không suy biến, giải thuật đơn hình trên đây sẽ hội tụ về phương 
án tối ưu sau một số hữu hạn lần lặp. 

Bằng sự thống kê người thấy rằng nói chung giải thuật đơn hình sẽ hội tụ với số lần 


lặp ít nhất phải là từ m đến 3m ( m là số ràng buộc ) . 


Giải thuật đơn hình cơ bản 


Xét bài toán quy hoạch tuyến tính chính tắc 
min/maxz(z) = cŸz 

Ax=b 

>0 


Giả sử rằng sau khi hoán vị các cột trong A ta chọn được ma trận cơ sở B thoả sự phân 
hoạch sau đây : 


A =[BN] 


T = [egcn]ư 


C 
œŸ = [ep#w]|ư 


Giải thuật đơn hình cơ bản được thực hiện như sau : 


a- Tính ma trận nghịch đảo B-1 
b- Tính các tham số : 

. Phương án cơ sở khả thi tốt hơn 
zp=B lb=b 


?.R:==2Ú 


rùgh 


. Giá trị hàm mục tiêu z() = cb#ø 
. Ma trận W= B-1N 
c- Xét dấu hiệu tối Ưu : 

ST `.) 
mu =óny CN” ẤL cCặy — cnN 
- Nếu Ếy < 0 thì kết thúc giải thuật với phương án tối Ưu là : 
#;B.'.= B'b=b 


#zx=0 


rùgh 


và giá trị hàm mục tiêu là : 


2P) ED4'p 

- Nếu tồn tại ẽ¿ € £y mà ẽ; > 0thì sang bước d. 

d- Xác định chỉ số của phần tử pivot trong ma trận 
. Xác định chỉ số cỘt s của pivot 

é¿ = max{e; >0 € 6y} 


Nếu Ñ;„ < 0 thì giải thuật dừng, bài toán không có phương án tối ưu. Ngược lại thì tiếp 
tục. 


. Xác định chỉ số dòng r của pivot 


min{ hi vŸ 0} =-È (ý = 1/2) Ÿn) 
Ns ÁN 


Phần tử Ñ;; trong ma trận được gọi là phần tử pivot 
Trong trường hợp bài toán min 

c- Xét dấu hiệu tối Ưu : 

Em = nỊn — cầBN = STỦ — caN 

- Nếu z1; >0 thì kết thúc giải thuật với phương án tối ưu là : 
#a =1: 1b=b 


#zx=0 


r¡ùgh 


và giá trị hàm mục tiêu là : 


2z(+) = càzø 


- Nếu tồn tại ; € £wy mà ẽ; < 0thì sang bước d. 

d- Xác định chỉ số của phần tử pivot trong ma trận ý 
. Xác định chỉ sỐ cỘt s của pivot 

c¿ = max{| ¿ | 6y < 0 € £w} 


Nếu Ñ;, < 0 thì giải thuật dừng, bài toán không có phương án tối ưu. Ngược lại thì tiếp 
tục. 


. Xác định chỉ số dòng r của pivot 


min{ N. xi 0} .a (=1. 70)) 


Phần tử N;; trong ma trận được gọi là phần tử pivot 
e- Thực hiện các hoán vị : 

. CỘt thỨ s trong ma trận N với cột thứỨ r trong ma trận B 
. Phần tử thứ s trong c?; với phần tử thứr trong c5 

. Biến xs trong z1 với biẾn xr trong #» 

f- Quay về (a) 


Ví dụ : Tìm phương án tối ưu cho bài toán quy hoạch tuyến tính chính tắc sau đây bằng 
giải thuật đơn hình cơ bản 


maxz(#) = 2x1 + #2 
#1— #a+#ạ7—= 3Š 
ø1 + 2xa + #4 = 6 

—#Z1 + 2xaạ + #zs = 2 

z; > 0(7 = 1,2,3,4,5) 


tí 


Ta có : 


t' È Ïi #4 di 3 
A=|l1 3# |010|b=|s 
ciP: .6: ` (1500910 71 2 
NB 
ØỞ=[zi 8a | 8s #84 #a| 
#N#p 
c7 =|2 1 | 00 0| 
CẶCp 
Lần lặp1 


a- Tính ma trận nghịch đảo B-1 


B1=B— 


© 
© CC 


0 
0 
1 
b- Tính các tham số 


. PhƯơng án cơ sở khả thi tốt hơn : 


%3 
+4 
%š5 


rùgh 


100 
010 
001 


rùgh 


+#bB— 


. Giá trị hàm mục tiêu : 
2#) = Con = |0:.0.:0)/1161:= 0 


. Tính ma trận : 


Lữ ñ[Ƒ3: ›=ĩ 1 ¿ii 
N=-EB`ÌN-|0 1 0l|1 2|-= 11 2 
027 JI<-f- =j, Ð 
c- Xét dấu hiệu tối Ưu : 
†: =ĩï 


ấy =cy—csWN=[2 1|-[0 0 0]|1 2|= 


Chuyển sang bước d 

d- Xác định chỉ số của pivot 

. Xác định chỉ sỐ cỘt pivot s : 

6; =max{e; > 0€ 6y} =max{2, I1} =2 = ei 


Vậy s=1 


ñÌ 
ị 
=l 


rùgh 
Ma trận cột s=1 trong ma trận Ñ là 


: b, H bị b; } h 3 6 bì 
m1n4 — = mln4 ——,—= = min4 ~-,—t—=ä37=_— 
{x.) nến, 1117 Nụ 


VẬậyr=1 

e- Hoán vị 

. Cột thỨ s=1 trong ma trận N và cột thứ r=1 trong ma trận B 
. Phần tử thứ s=1 trong c1; với phần tử thứ r=1 trong c? 


. Biến thứ s=1 trong z„ với biến thỨ r=1 trong z» 


=3. f1 30 .0 RseL:Ì 2E: f0 
| tr. Ø.10t04140|2/6r|0-36/ 7l T40 
-1 9 |001 0 |—1 01 
chh= T5 9Ÿ ^t5Pi0)926/cFee T051 f8) 
[ri đa. | #3 đụ #s] >XÍ =[#a4 #øạ | #1 #4 | 


f- Quay về bước a 
Lần lặp 2 


a. Tính ma trận nghịch đảo B-1 


b- Tính các tham sỐ 
. PhƯơng án cơ sở khả thi tốt hơn : 


+1 
4 
%š5 


righ 


righ 


. Giá trị hàm mục tiêu : 


3 
z(z) =cbzp =[|2 0 0]|3| =6 
5 


. Tính ma trận : 


N=B ÌN- |1 1 0|Ì0 2| =|-1 3 


c- Xét dấu hiệu tối Ưu : 


II 
Zy =cyT— cgsÑ =[0 1J—[200]|-1 3 | =[-23| 
Ï.. SĨ 


Chuyển sang bước d 
d- Xác định chỉ số của pivot 
. Xác định chỉ sỐ cỘt pivot s : 


c; = max{ey > 0 € £„y} = max{3} = 3 = ca 


Vậy s=2 
-1 
b) 
rùgh 
Ma trận cột s=2 trong ma trận Ñ là 
[l 
Na—= 


. Xác định chỉ số dòng pivot r : 


: b; : bạ by xo b; 
min4 — — min4 ——,— = min4 ,—‡ —=l=—= 
{z } {các} {3:1} N; 


1S 


Vậyr=2 
e- Hoán vị 


. Cột thứ s=2 trong ma trận N và cột thứ r=2 trong ma trận B 


. Phần tử thứ s=2 trong c1 với phần tử thứ r=2 trong c? 


. Biến thứ s=2 trong z1 với biến thứ r=2 trong + 


1 1: | 1⁄0. J. St: s10 
A=l|0 2 | 1 10 >A=l01| 1 2 0 
ñ- 3 lÌ: SE 3Ô ñ' Ú T61. 7 đ 


cT=|0 1 |2 0 0+cT=l0 0 |2 1 0| 

#1 = [zs %3 | #1 #4 mg] —> x” = [#3 #x.° || -3#t. '#ð° si 
f- Quay về bước a 

Lần lặp 3 


a. Tính ma trận nghịch đảo B-1 


b- Tính các tham sỐ 
. PhƯơng án cơ sở khả thi tốt hơn : 


#1 
+2 
%5 


r2nh, 


“. .9”'eẻ 


r+gh 


r+gh 


r+gh 


+ã 


#4 


r+gh 


r¡ùgh 


. Giá trị hàm mục tiêu : 


4 
Z1) <= G 6n =Ï2 1/0) =0 
4 


. Tính ma trận : 


lÌÌ 
N=B'1N= 


c- Xét dấu hiệu tối Ưu : 


ẩyạ„=œy—c„N=|0 0|— [210] 
2 
3 


: dừng 


llÌ 


Vậy phƯơng án tối ưu sẽ là : 


Giá trị hàm mục tiêu là z(x) = 9 với x1 = 4 và x2 = 1 


Chú ý trong trường hỢp suy biến 


Trong trường hợp bài toán suy biến, nghĩa là b„ = 0, ta có : 


cho nên giá trị của hàm mục tiêu không thay đổi khi thay đổi cơ sở, vì : 
2(“)].=:2®)-<.baa= 20) 

Vậy thì, có thể sau một số lần thay đổi cơ sở lại quay trở về cơ sở đã gặp và lặp như 
vậy một cách vô hạn. Người ta có nhiều cách để khắc phục hiện tượng này bằng cách 


xáo trộn một chút các dữ liệu của bài toán, sử dụng thủ tục từ vựng, quy tắc chọn pivot 
để tránh bị khử: 


GIẢI THUẬT ĐƠN HÌNH CẢI TIẾN 
Một cách tính ma trận nghịch đảo 
Trong giải thuật đơn hình cơ bản hai ma trận kề B và chỉ khác nhau một cột vì vậy có 


© ^ . 2 1 ^ > b$ A; ` Z1 . z 2 ^ ^ 
thể tính ma trận nghịch đảo Ö_ một cách dễ dàng từ B-1 . Để làm điều đó chỉ cân nhân 
(bên trái) B-1 với một ma trận đổi cơ sở được xác định như sau : 


—¬ 
©= 

= 
©= 


Ụ 1 = 0 
SP (Út đấy =— . 0| >dòngr 
0 0 — lj 
† côtr 
Khi đó : 
~1=„B1 
B 


Ta thấy rằng ma trận đổi cơsở được thiết lập giống như một ma trận đơn vị mxm, 
trong đó cột r có các thành phần được xác định nhƯ sau : 


=>: đối với thành phẩn¡ r. 
s-: đối với thành phần r. 
Khi mà ma trận cở sở xuất phát là ma trận đơn vị, sau một số bước đổi cơ sở B0 B1 B2 


NgàẫC Bq tương ứng với các ma trận đổi cơ sở 0 1 2....... q-1 người ta có cách tính ma 
trận nghịch đảo như sau : 


Quy hoạch tuyến tính dạng chuẩn 


Quy hoạch tuyến tính dạng chuẩn là quy hoạch tuyến tính chính tắc mà trong đó có thể 
rút ra một ma trận cơ sở là ma trận đơn vị. Quy hoạch tuyến tính chuẩn có dạng : 


min/maxz(z) = cŸz 


[JN]+ = b 
zẩ>0 


Giải thuật đơn hình cải tiến 


Từ những kết quả trên người ta xây dựng giải thuật đơn hình cải tiến đối với bài toán 
qui hoạch tuyến tính (max) dạng chuẩn nhƯ sau : 


a- Khởi tạo 
Ao=4A 
bạ =b 


b- Thực hiện bước lặp với k = 0,1,2, ... 


. Xác định phƯơng án cơ sở khả thi : 


. Tính giá trị hàm mục tiêu : 


z(z*) = c?,#p, = ch, 


. Xét dấu hiệu tối Ưu : 


+YB, — Dư 
ỞŸN, — 0 
rùgh 


là phương án tối ưu 


z(z*) = Có ng = cp,b; là giá trị hàm mục tiêu 
- Ngược lại thì sang bước (c) 

c- Cập nhật các giá trị mới : 

.Tính pivot 

.Tính ma trận chuyển cơ sở k 

.Tính Akai = HhẬt 

.Tính bg¿ị¡ = *b; 
.Tăng số lần lặp k=k+1. 
Quay về bước b 

Ví dụ 


Giải bài toán quy hoạch tuyến tính sau đây bằng phương pháp đơn hình cải tiến : 


max Z(#) = 2xị + #2 
#1 — #2 T†#a3a =3 
#1  2xạ #4 =Ö 
—Z1 + 2xa + #øg = 2 


ø; > 0(j = 1,2,3,4,5) 


t{ 
Bước khởi tạo 
=ƒ /l*$b 028 3 
PU 2 | 01 0lỗạ= |6 
S1: đốc jllr 8-2091 2 
NoBọ 


Bước lặp k=0 


+3 
%4 
%5 


=0 
#Ng 


righ 


= |00 0] 


= 21000] — |0 00] 
1-1100 
12010 
—12001 


rùgh 


llÌ 


ST 0 092 TẢ T5 
cọ =c” —ch,Áo 


Ï 
—] 
righ 
| 
3 Suy ra pivot : 1i = 1 
6 
P, 
r¡igh 
[ll 
L 
—= le 1 
1 


11100 1-1100 
12010 03-110 
—12001 0O1I10I 


rìgh righ 


[Ì [Ì 


1 0 0|r¿gh r¡ígh 
bạ = u°bg =|—-l I1 0 = 
XI G lân: 


II  H 


Bước lặp k=1 


tẾW =='Ú 


righ 


= |20 0] 


1-1100 
03-110 
01101 
rìgh 
ếi =c” —cp 4i =[21000] — [200] 


llÌ 


=[03-200] 


áy 250) SH DÌNHDS ĐiÁe 9 
3 
bì 
rìgh 
[Ï 
¬ 
ứ=|0 ‡ 
0 " 1 


1:ili0n8c v0.0 
03-110 01-41 


3 3 
Ø1. 1,01 4 1 
00-1 


l II 


Š 
5 4 
: : 0 righ righ 
0 —# ïÌ 


II djI 


righ 


SN ==Ú 


righ 


= [210 


sẽ r¡igh 
đ; =cT —cp? 4; = [21000] — [210] 


[Ì 


=[00 -1 -10]: thoả dấu hiệu tối ưu. 


Vậy kết quả của bài toán là : 


“m> €Œ CC 


. Phương án tối Ưu x = x2 = righ 


[II 


. Giá trị hàm mục tiêu z(x) = 9 


Phép tính trên dòng - Bẳng đơn hình 

Các bước thực hiện giải thuật đơn hình cải tiến được trình bày lần lượt trong các bảng, 
gọi là bảng đơn hình. Trong thực hành, để cập nhật những giá trị mới ta có thể làm như 
sau : 

. Tìm pIivot. 

. Chia dòng chứa pivot cho pivot. 

. Khử các phần tử trên cột chứa pivot. 

. Tính dấu hiệu tối ưu. 


. Tính giá trị hàm mục tiêu . 


CBụ ?Bụ #1 #2 #3 %4 %5 bọ 


CBị 


CB; 


?Bì 


2B, 


+1 


Z1 


+2 


Z2 


+3 


-2 


Z4 


%5 


%5 


z(x0) 


Z(x1) 


Z(x2) 


PHƯƠNG PHÁP BIÊN GIÁ CẢI BIÊN 

Bài toán cải biên 

a- Cải biên bài toán quy hoạch tuyến tính 

Người ta có thể biến đổi một bài toán quy hoạch tuyến tính chính tắc thành dạng chuẩn 
bằng cách cộng một cách phù hợp vào vế trái của ràng buộc ¡ một biến giả xnti 0 để 
làm xuất hiện ma trận đơn vị. Vì các biến giả cải biên có ảnh hưởng đến hàm mục tiêu 
nên cũng sẽ có sự cải biên hàm mục tiêu. 


Vậy, người ta có thể biến đổi bài toán quy hoạch tuyến tính tổng quát, gọi là bài toán 
xuất phát, thành bài toán dạng chuẩn, gọi là bài toán cải biên (mở rộng) 


Ví dụ : 
Biến đổi bài toán quy hoạch tuyến tính sau đây thành dạng chuẩn 


maxz(z) = 2XxỊI †#a + #3 — Z4 
ø1 + 5xa + 5xa = 25 
—4xa — #a + Öxa = l8 
3x›¿ + 8xa = 28 


z¡ > 0( = 1,3,3,4) 
tí 


Bài toán xuất phát có các biến, ma trận ràng buộc và chỉ phí : 


1505 
0-4-16 
0306 


r¡ùgh 


An: = [#ixaxzxa] [ [ 


A = 


Bằng cách thêm biến giả x5, x6 lần lượt vào ràng buộc 2 và 3. Ta được bài toán cải 
biên : 


max{z () = 2xị + #s + # — #a — M(œs + #¿) 
#1 + 5xạ + 5xa = 25 
—4xa — #3 + xa + zs = 18 
JdXa + Öxa + øạ —= 28 


z¡ > 0(7 = 1,2,3,4,5,6) 
t 


z (œ) là hàm mục tiêu cải biên sẽ được giải thích trong phần tiếp theo. 


Các biến, ma trận ràng buộc các hệ số và chi phí của bài toán cải biên là 


150500 
0-4-1610 
030601 


righ 


z”. = [#1XaXzXaX;Xö] [ [ 


A = 


b- Quan hệ giữa bài toán xuất phát và bài toán cải biên 
Người ta kiểm chứng rằng : 


- Nếu #z# = [z#a...z„| là phương án (tối ưu) của bài toán xuất phát thì 
#T = |ztza...z„0 0...0] là phương án (tối ưu) của bài toán cải biên tương ứng. 


Vậy nếu bài toán cải biên không có phương án tối ưu thì bài toán xuất phát cũng sẽ 
không có phương án tôi ưu. 


- Nếu #? = [z1za...z„0 0...0] là phương án tối ưu của bài toán cải biên thì 


œ? = |#zs...„] là phương án tối ưu của bài toán xuất phát 


- Nếu bài toán cải biên có một phương án tối ưu mà trong đó có ít nhất một biến giả có 
giá trị dương thì bài toán xuât phát không có phương án tối ưu. 


- Nếu bài toán cải biên (dạng chuẩn) có phƯơng án tối ưu thì cũng sẽ phương án cơ sở 
tối ưu. 


Ví dụ 


1- Xét bài toán : 


minz(#) = #1 + 2xa + #4 — ðXg 
—3xs — xạ = 0 
#2 — Xa — Xa T— 2X; — 5 


1 2 4 1 ¬- 
#1i— g2 † 33T #4 7s =$ 


ø¡ > 0(7 = 1,2,3,4, 5) 
{ 


Bài toán cải biên không có phƯƠơng án tối ưu nên bài toán xuất phát cũng không có 
phương án tối ưu. 


2- Xét bài toán : 


mỉn z() —= —lồxị + 7xa + 9x; 
—ÝZ1 = 32 + #g —= 3 
—ÖXỊ + he) =ĩ 


{ 


Phương án tối ưu của bài toán cải biên : 
kết: ⁄#zc øs: 24]: |0 H { 0| 
Phương án tối ưu của bài toán xuất phát : 


[#1 Z2 #s] = |0 H TH 


3- Xét bài toán : 


mỉn z() = 2xị + 4xa — 2x; 
ø1 — 2Xa + #a = 27 
2xi +#a + 2xs = 50 
#1 — 8a — #øs < 18 


#j( = 1,2,3) 
t 


Phương án tối ưu của bài toán cải biên : 
ly: 3. 37 28x: ' 2867 406]. = |0; 0 35.:43. 3: Dị 
Bài toán xuất phát không có phương án tối ưu. 


Hai phƯơng pháp biến giả cải biên thương dùng là phương pháp hai pha và phương pháp 
M vô cùng lớn . 


Phương pháp hai pha 

Pha 1 

Tìm phương án tối ưu cho bài toán cải biên với hàm mục tiêu cải biên là : 

min (tổng tất cả biến giả cải biên) 

Pha 2 

Tìm phương án tối Ưu cho bài toán xuất phát với phương án cơ sở khả thi xuất phát là 
phương án tối ưu tìm được ở pha 1. Ở pha 2 này các biến giả cải biên bị loại ra khỏi ma 
trận các hệ số ràng buộc, và vectơ chi phí được cập nhật lại, do đó dấu hiệu tối ưu cũng 
được cập nhật lại 

Đây là phương pháp thuận lợi cho việc lập trình ứng dụng giải thuật đơn hình cải tiến. 


Ví dụ : Xét bài toán quy hoạch tuyến tính 


maxz(#) = 3x + 4xa + Z3 
#1 + 2xạ + 2xs < 


©c2|—1 c2|œ 


#1 + 2x; + 3x; > 


{ 


Đưa bài toán về dạng chính tắc bằng cách thêm biến phụ x4, x5 ta được 


max2z(ø) = 3xi + 4xa + #¿ 
#1 + 2xa + 2xa + Z4 — 


©cz|—4 ©|œ 


ÿ1-E2Xãz 11 3Xã —ứnp= 


z;¡ > 0(7 = 1,2,3,4,5) 
{ 


Ma trận các hệ số ràng buộc là : 


12210 
1230 —1 
righ 
A= không chứa ma trận đơn vị 


Áp dụng phương pháp đơn hình cải biên hai pha nhƯ sau : 


Pha 1 


Thêm biến giả (cải biên ) x6 0 vào ràng buộc thứ hai để được ma trận đơn vị . Khi đó 
bài toán cải biên có dạng : 
minu(#) = #s 
#1 + 2x2 + 2x3 + 4 = Š 
#1 + 2xa + ŠXa — #gs + ø — 


©|~1 


z¡ > 0(7 = 1,2,3,4,5,6) 
{ 


Có ma trận các ràng buộc là : 


122100 
1230 — 11 
rùgh 
có chứa ma trận đơn vị 
| 
A-—= 


Giải bài toán cải biên bằng giải thuật đơn hình cải tiến 
Khởi tạo 


— D227202si; “JLC 2U 
Ao= bọ = 
1/292: 509- J ỦY, “SEN s 


©œ|—t c2|œ 
=“.=—— 


c®=[0 0000 1] 


Bước lặp k=0 


CBụ ?Bạ x1 x2 x3 x4 x5 


0 4 1 2 5) 1 0 
1 6 1 2 3 0 -1 
cT 0 0 0 0 0 
Gà s† -2 -3 0 1 


CB; ?Bì x1 x2 x3 x4 X5 
0 4 _ Ỹ 0 1 ? 
0 3 3 Ỹ ï 0 —# 
cT 0 0 0 0 0 
Eh 0 0 0 0 0 


Ta được phƯơng án tối ưu . Xong pha 1 . Chuyển sang pha 2. 
Pha 2 
Loại bỏ biến giả cải biên x6 0 


Khởi tạo 


©32|2 c2|b2 
œœ 
| ©|2 
c|— 

LE—====ễ—=-i 


‹©|~t ©ÍbB _ C© 


ai 
= 
| 
EE—c=—=ee.mđ 
_——c¬ 


cT=[3 4 1 0 0| 


x6 


w(0) 


Bước lặp k=0 


Cụ 1? Đụ 
0 4 

1 3 
cT 

Số 


CD; ? Dị 
0 4 
4 2 
cT 

cï 


Bước lặp k=2 


CB, 2B, 


x1 


có|— 


x1 


}.ớ2) x3 x4 
Ỷ 0 1 
$ 1 0 
4 1 0 
bi 0 0 
x2 x3 x4 
0 -1 1 

J Ỷ 0 
4 † 0 
0 -Ð 0 

x2 x3 x4 


x5 


©2|k2 


¬ 


| 
ko|— 


x5 


Bước lặp k=3 


CB; 


1B; 


0 0 -1 
1 
3 1 1 
&) 4 1 
1 0 -3 
x1 x2 x3 
0 0 -1 
1 2 2 
&) 4 1 
0 -2 -D 


Kết quả của bài toán đã cho : 


. PhƯơng án tối ưu 


đi =Ỹ 
1= Ú 
Xu =0 
Ứ4:= 
#g = $ 
tít 


. Giá trị hàm mục tiêu z(x)=z(x3)= 8 


Phương pháp M vô cùng lớn 


x4 


x5 


Phương pháp M vô cùng lớn ( M là số vô cùng lớn ) tương tự như phương pháp hai pha, 
ngoại trừ Ở pha 1 hàm mục tiêu cải biên có dạng sau đây cho bài toán max/min 


max [z(x) - M*( tổng các biến giả cải biên) ] 
min [z(x) + M*( tổng các biến giả cải biên) ] 


Bằng phương pháp này, trong quá. trình tối ưu, các biến giả cải biên sẽ được loại dần ra 
khỏi ma trận cơ sở : tất cả đều bằng 0. Nếu trong quá. trình tìm phƯƠơng án tối ưu mà 
không loại bỏ được các biến giả cải biên ra khỏi cơ sở thì bài toán vô nghiệm. 


So với phương pháp hai pha thì phương pháp này tránh được việc phải cập nhật lại dữ 
liệu cho bài toán gốc nhưng không tiện lợi bằng trong lập trình ứng dụng. 


Ví dụ : Xét bài toán tương tự như trên 


maxz(#) = 3xị + 4xạ + Z3 
#1 + 2xa + 2xa + #4 — 


#1-+2xa-+ 3Xa — #y = 


c2|—1 c2|œ 


z¡ > 0(7 = 1,2,3,4,5) 
{ 


Thêm biến giả cải biên x6 0 vào ràng buộc thứ hai đồng thời cải biên hàm mục tiêu 
theo như trên ta đƯỢc : 


maxu(z) = äxị + 4xa +zs — Mxs 
#1 + 2xa + 2x + z4 —= Ỹ 


#1 + 2x2 + 3X; — Øg + #a — {+ 


z¡ > 0(7 = 1,2,3,4,5,6) 
{ 


Tìm phương án tối Ưu cho bài toán cải biên này bằng phương pháp đơn hình cải tiến 


Khởi tạo 


" 1 2 2 1 ` 
TT tJ5=| 


c=[3 4 1 0 0 —MỊ 
Í.2)-Số, (0<: 


©›|—t ©2|@œ 


Bước lặp k=0 


CBụ tụ #1 #2 %3 #4 #5 %6 bọ 
0 4 1 2 2 ] 0 0 š 
M6 1 2 3 0 -1 1 3 
cT 3 4 1 0 0 -M w(0) 
sã, 3~_M 422M 113M. 0 M0 _= 


CB, 12B, #1 #2 %3 %4 #5 %6 bị 
0 4 75 | 0. |7 |v —' H 
c3 | LỆ 1 9  -‡ LỆ h 
cT & 4 1 0 0 -M w(x1) 
aT h { |0 0 _ —.M 5 


Do x6 = 0 (vì ngoài cơ sở) nên bị loại ra khỏi bảng và ta tiếp tục tìm phương án tối ưu 
cho bài toán gốc đã cho có phƯơng án cơ sở khả thi được khởi tạo như sau : 


CBụ ?Bì Z1 Z2 #3 4 #5 bọ 


' 4 $ |‡ |0 Ị : # 
Ị 3 H Ỹ 1 0 -‡ : 
cT 3 4 ƒ 0 0 z(x0) 
ấị ‡ÿ |J# |° |° |3 Ỹ 


Các bước tiếp theo được thực hiện giống như phương pháp hai pha. 


QUY HOẠCH TUYẾN TÍNH SUY BIẾN 


Khi thực hiện thuật toán đơn hình trường hợp bất thường có thể xảy ra là khi xác định 
biến ra thì tồn tại tỷ số h = 0, tức là tồn tại bi=0, hay không có tỷ số nào dương thật 
sự. Người ta xem đây là trường hợp suy biến. Khi một bảng đơn hình rơi vào tình trạng 
suy biến thì có thể gây khó khăn mà cũng có thể không khi ta tiếp tục thực hiện thuật 
toán đơn hình. 


Các ví dụ về quy hoạch tuyến tính suy biến 
Ví dụ 1 : xét quy hoạch tuyến tính : 


mỉnz(#) = 7 + # — #2 
#1 — 2xạ <S 2 
—jxi <Ổ 
—2xị <0 


%1;2 > 0 


tÍ 


Đưa bài toán về dạng chuẩn : 


mỉnz(#) = 7 + #1 — #2 


#1 — 2Xa + za = 2 


—ðXi + #4 7= 6 
—2XỊ + #s—= 0 
#1,#a > 0 
lại 


với ma trận hệ số là : 


x1 Xã x3 x4 x5 
lj -2 1 0 0 
-3 0 0 1 0 
-2 0 0 0 } 


có chứa ma trận đơn vị. Áp dụng thuật toán đơn hình cải tiến ta được : 


cB 1B x1 x2 xó x4 x5 
0 3 1 -2 1 0 0 
0 4 -3 0 0 lị 0 
0 s -2 0 0 0 1 


W=7 


Đây là trường hợp suy biến, biến vào là x2, nó được tăng lên đến mức vẫn thỏa những 
điều kiện về dấu của các biến trong cơ sở x3, x3, x5. Đó là : 

zsza=2+2xạ>0 #z>—Ÿ 

#a=6+0xa>0 Vz¿ạ >0 

#sa =0 +0xạ>0 Vza>0 


tÍ tÍ 


Như vậy x2 có thể lớn tùy ý nên hàm mục tiêu không bị giới nội. Vậy bài toán không có 
phương án tối ưu. Trường hợp này ở bảng đơn hình không có tỷ số nào dương thật sự 
để xác định biẾn ra. 


Ví dụ 2 : xét quy hoạch tuyến tính : 


mỉnz(#) = 7 + #1 — #2 
#1 + 2xạ < 2 
—jxi <Ổ 
—2XỊ < 0 


fj405:-” 


lht 


Đưa bài toán về dạng chuẩn : 


minz(#) = 7 + #1 — #2 


#1 + 2xạ + #s = 2 


—ðXi + #4 —= 6 
—2XỊ + #s—= 0 
#1,#a > 0 
lại 


với ma trận hệ số là : 


x1 Xã 
lj 2 
-3 0 
-2 0 


có chứa ma trận đơn vị. 


cB 1B 
0 ° 
0 4 
0 Đ 


x3 x4 x5 
1 0 0 
0 1 0 
0 0 1 


Áp dụng thuật toán đơn hình cải tiến ta được : 


x1 x2 xó x4 x5 
1 2 1 0 0 
-3 0 0 lị 0 
-2 0 0 0 lÌ 


W=7 


cT 1 sỉ 
eT ï vÍ 
cB 1B x1 
1 
si Đ) 3 
0 4 -3 
0 5 =>, 
cT 1 vĩ 
=T 3 
Ỹ 0 


Đây là bảng đơn hình tối ưu. 


Ví dụ 3 : xét quy hoạch tuyến tính : 


x2 


k2|D— 


x4 


x5 


minw(z) = -ä3 + Z1 — 2xa + 3g 
Z1 + +3 <1 


—#+ +#a—#s<0 


#1,#a2,#s > 0 


{ 


Đưa bài toán về dạng chuẩn : 


mỉnw(#) = -3 + 2ø — 2x; + 3 
1 + 33 + 4 =i 


—#1 + #a — #s + z; =0 


#19021; w„0s >>.U 


{ 


với ma trận hệ số : 


x1 x2 x3 x4 x5 
1 1 

5 0 5 1 0 

-1 1 1 0 1 


có chứa ma trận đơn vị . Áp dụng giải thuật đơn hình cải tiến : 


cB 1B 

0 4 

0 5 
cT 3 
e7 3 


X2 vào, x5 ra 


cB 1B 

0 4 

-2 2 
XE 
2T = 


x1 vào, x4 ra 


2|© 


.) 


x2 


2|©2 


t2|© 


x2 


t2|© 


k2|— 


k2|— 


x4 


x4 


0 

1 

w=-3 
x5 


w=-3 


cB 1B x1 x2 x3 x4 

3 1 1 0 1 2 

-2 2 0 1 0 2 
cr 3 -2 Ỷ 0 0 
TH 0 0 1 3 Ũ 


Đây là bảng đơn hình tối ưu 
Ví dụ 4: xét quy hoạch tuyến tính 


minz(z) = —10Z + 57z¿ + 9x + 24z¿ 
U,5#1 — 5,5; — 2,Bax -+ 9x < Ö 
0,571 — 1,5#a — Ö,ðzs + za¿ < 0 
#øiSl 


11;#2;#Z3;Z4 > 0 


tÍ 


Đưa bài toán về dạng chuẩn 


-6 


minw(#) = —10zq + ð7za + 9xạ + 24z4 
0,521 — 5,52 — 2,523 + 9x4 + #s — 0 
0,5Z1 — 1,5#a — Ö,ðZ3ä + #4 + #6 = 0 
#1 +; = Ì 


%1;#2;13;4;5;6;#7 > 0 


t{ 
với ma trận hệ số 
x1 x2 x3 x4 x5 x6 X7 
0,5 -Ð,D -2,5 9 1 0 0 
0,5 -1,5 -0,5 1 0 1 0 
1 0 0 0 0 0 1 


có chứa ma trận đơn vị . Áp dụng phương pháp đơn hình cải tiến 


cB 1B x1 x2 x3 x4 x5 x6 
0 5 0,5 -D,D -2,5 9 1 0 
0 6 0,5 -1,5 -0,5 1 0 1 


cT -10 57 9 24 0 0 0 


eT -10 57 9 24 0 0 0 


x1 vào, x5 ra 


cB 1B x1 x2 x3 x4 x5 x6 
-10 1 Ỳ TIM -5 18 2 0 
0 6 0 4 2 -8 Ti 1 
0 ớ/ 0 11 5 -18 -2 0 
c. 10 57 9 24 0 0 0 
G* 0 -53 -41 204 20 10 0 


X2 vào, x6 ra 


cB 1B x1 x2 x3 x4 x5 x6 
-10 1 1 0 0,5 -4 -0,75 2,75 


S7 ¿ 0 J 0,5 -2 -0,25 0,25 


X3 vào, xÍ ra 


cB 


S7 


1B x1 

3 2 

2 -1 

7 1 
-10 DĐ 
29 0 


X4 vào, x2 ra 


x2 


24 


98 


x3 


0 


6,75 


x4 


0 


13,25 


-2,75 


1 


X7 


X5 vào, X3 ra 


cB 


24 


1B 


x3 


x4 


x4 


70,5 


x5 


x6 


X7 


X7 


x6 vào, x4 ra 


cB 1B xi xa x3 x4 x5 x6 Xơ b 
0 5 0,5 -5,5 j5 9 I 0 0 0 
0 6 0,5 ST -0,5 1 0 ƒ 0 0 
0 7 1 0 0 0 0 0 1 1 
w=0 
c -10 57 9 24 0 0 0 
G -10 57 9 24 0 0 0 


Bảng đơn hình hiện thời giống với bảng đơn hình xuất phát : đây là hiện tượng xoay 
VÒng . 


XỬ lý trường hợp suy biến 


Theo các ví dụ trên, trong trường hợp quy hoạch tuyến tính suy biến thì sau một số lần 
lặp có thể phương án nhận được vẫn như cũ mà không có sự thay đổi nào, có thể 
phương án nhận được tốt hơn, có thể phương án nhận được là một phương án đã nhận 
trước đó rồi và từ đó cứ xoay vòng mãi. Do đó nếu không có biện pháp phòng ngừa thì 
thuật toán đơn hình sẽ có thể kéo dài vô tận. 


Khi thực hiện thuật toán đơn hình thì hiện tượng suy biến xảy ra khi có sự tình cờ khử 
lẫn nhau làm cho tồn tại b; nào đó bằng 0. Trong trường hợp này có thể có nhiều biến 
thỏa điều kiện của biến ra. Gặp trường hợp này cần phải lựa chọn biến ra sao cho tránh 
được hiện tượng xoay vòng. 


Người ta thường dùng phương pháp nhiễu loạn, phương pháp từ vựng để tránh sự tình 
cờ khử lẫn nhau này. Trong thực tiển tính toán người ta đã đề ra một quy tắc xử lý khá 
đơn giản, gọi là quy tắc Bland, khi dùng giải thuật đơn hình giải các quy hoạch tuyến 

tính suy biến, đó là : 


Với xk là biến vào , biến ra xr được chọn là biến có chỉ số nhỏ nhất thỏa điều kiện chọn 
biến ra : 


min{ bị ;ữïk > 0} = 1n) 


địk 
Ví dụ : 
Xét quy hoạch tuyến tính suy biến : 


minw(#) = —34 — +2*xg — #s + lỗ#; 
#\ + 34 — 2Xg — øe + 12; = Ö 
#2 + 24 — #s — $e + 2ø; = 0 


#3 + #s + øg — 9x; = 2 


%1;2;#3;⁄4;5;Z6;#7 c. 0 


tí 


Áp dụng quy tắc Bland ta thấy : 


cB 1B x1 x2 x3 x4 x5 x6 x7 b 
0 1 1 0 0 3 -2 +1 12 0 
0 2 0 j 0 D +1 na $Ỹ 0 
0 3 0 0 1 0 1 1 -9 2 


Biến ra có thể là x1 hay x2 


cB iB 
A 

mẽ 4 

0 3 

0 3 
cT 0 
—T 4 

Biến ra là x2 
cB 1B 


| 
2|© 


_4 
3 
. Chọn x1 
x2 x3 
0 0 
1 0 
0 1 
A 
>Sg) 
0 
x2 x3 
&) 0 


x4 


x5 


-6 


x5 


16 
16 
x6 X7 
-3 36 
4 34 
g >— nay 
1 -0 
w=0 
16 
64 
x6 X7 
1 2 


| 
©2| 


Biến ra có thể là x4 hay x5 . Chọn x4 


cB 1B x1 


-1 6 —Š 

2 5 + 

0 3 n 
cT 0 0 
eT -2 6 


Biến ra là x5 


x3 


| 
©2| 


— 


x5 


©|k2 


c|— 


16 


30 


x6 


16 


32 


w=0 


X7 


w=0 


cB 


cT 


1B 


Biến ra là x3 


cB 


cT 


1B 


x1 x2 
0 -6 
1 6 
0 6 
0 0 
-6 0 
x1 x2 
54 
: 31 
18 
1 T 
6 
0 3T 
0 0 


x3 


©|t> 


x5 


x6 X7 
1 -40 
0 -28 
0 J1 
W_—= 
16 
-24 
x6 X7 
1 0 
0 0 
0 1 
48 
— 31 
16 


eT “2 24 — 187 128 
Đ 31 31 93 31 0 0 


Đến đây không còn hiện tượng suy biến. 

Biến vào là x7 

CÂU HỎI CHƯƠNG 2 

1- Trình bày cơ sở lý thuyết của thuật toán đơn hình cơ bản. 

2- Định nghĩa quy hoạch tuyến chuẩn. 

3- Trình bày các bước lập bảng đơn hình theo phép toán trên dòng . 


4- Cải biên một quy hoạch tuyến tính tổng quát như thế nào ?. Cách giải quy hoạch 
tuyến tính cải biên và quy hoạch tuyến tính gốc. 


BÀI TẬP CHƯƠNG 2 
1- Tìm phương án tối Ưu của bài toán sau đây bằng phương pháp đơn hình cơ bản 


m1nz = -2xỊ — 2Xxa 
maxz = 3x + 2x2 


2xi +za<4 
-Xi+#a< 4 ' ' 

2XỊ + 3X2 < h) 
#1 + 2xa < 14 

4xi +#a <5 
ðxị + 2x¿ < 30 

#1 +ðxa <1 

01;ø>.10 
a)- b)- Z1,zz¿ >0 


tí (tt 


mìn w = #¡ + 2% + #s 
#1 + #2 + #3 + #4 + #s = ð 
#2 + #3 + #4 — #s —= 2 
#3 — øa + øg = Ì 
#1;Xa,Xs,Xa,Xg > Ö 


c)- 


tí 


2- Tìm phƯơng án tối ưu của bài toán sau bằng phương pháp đơn hình cải tiến 
a)max z = 5x1 + 3x2 

2x1I+2⁄2 90 

XI:‹.U 

x1,x2 0 

b)max z = x1 + 2x2 

2X PoOX2. ý 

x1I-x2 1 

x1 0,x2 0 

c) max z = 5X1 + 3x2 + x3 

2x1+3x2-x3 4 

3xI-x2+2x3 2 

XI} x2Z huïa 3 

x1I 0,x2 0,x3 0 

3- Tìm phƯơng án tối ưu của các bài toán sau bằng phương pháp biến giả cải biên. 
a)max z = 3x1 - x2 


2x1+x2 100 


x1 10 


b)min w = 3x1 + x2 
xXI'tx dã 

2Xl' 9 

x1,x2 0 

c)max z = 3x1 + x2 - 3x3 


xi+2x2-x3=2 


-10x2 + 5x3=5 
-3Xx2+2x3=4 
xi 0, I=13 


mInW = -XỊ — 3Xa 
maxz = 2xỊ -+ Ôxa 


#i+za>3 
—#] — #2 — #ạ S —2 
=1 đa = 
2XỊ — #2 + Z3 < Í 
#1 + 2xạ < 4 
Ø1,#2,#s > 0 
#1;2. >U 


d)- e)- 


tụ ( 


maxz = z¡ +3xa „minw = 2xỊ +#a 


—#1 — #¿ S —3 —#1 + #¿ <€ —1 

—=#++#a S—l =2 S2 

—# + 2xa < 2 si 
#1,#a > 0 #1;#a > 0 


Ð- 8)- 


tí tí 


Bài toán đối ngẫu 

Chương này trình bày trình bày khái niệm đối ngẫu, các quy tắc đối ngẫu và giải thuật 
đối ngẫu. Đây là các kiến thức có giá trị trong ứng dụng vì nhờ đó có thể giải một quy 
hoạch tuyến tính tỪ quy hoạch tuyến tính đối ngẫu của nó. 

KHÁI NIỆM VỀ ĐỐI NGẪU 

Đối ngẫu là một khái niệm cơ bản của việc giải bài toán quy hoạch tuyến tính vì lý 
thuyết đối ngẫu dẫn đến một kết quả có tầm quan trọng về mặt lý thuyết và cả mặt 
thực hành. 

Đối ngẫu cỦa quy hoạch tuyến tính dạng chính tắc 

Xét một bài toán quy hoạch tuyến tính dạng chính tắc 

về 


th 2(2) SE 


Giả sử rằng x* là phương án tối ưu cần tìm của bài toán và x0 là một phương án của bài 
toán thì mỘt cận trên của giá trị mỤc tiêu tối ưu được xác định vì : 


clx*£cTIx0 

Tuy chưa tìm được phương án tối ưu x* nhưng nếu biết thêm được một cận dưới của giá 
trị mục tiêu tối ưu thì ta đã giới hạn được phần nào giá trị mục tiêu tối ưu. Người ta ước 
lượng cận dưới này theo cách như sau : 

Với mỗi vectơ xT = [x1 x2... xn] 3 0 thuộc Rn chưa thoả ràng buộc của bài toán, tức là 
b-Ax!0 

người ta nới lỏng bài toán trên thành bài toán nới lỏng : 

min L(x,y) = cTx + yT(b - Ax) 

x=JU 


yT =[ y1 y2... ym] tuỳ ý Ì Rm 


Gọi g(y) là giá trị mỤc tiêu tối ưu của bài toán nới lỏng, ta có : 
g(y)= min { cTx + yÏI(b - Ax) } (x30) 
£ cIx + yÏI(b - Ax) 


Trong trường hợp x là phương án của bài toán ban đầu, tức là : 


b-Ax=0 
thì 
g(y) £ cTx 


Vậy g(y) là một cận dưới của giá trị mục tiêu bất kỳ nên cũng là cận dưới của giá trị mục 
tiêu tối ưu. 


Một cách tự nhiên là người ta quan tâm đến bài toán tìm cận dưới lớn nhất, đó là : 

max g(y) 

y tuỳ ý Ì Rm 

Bài toán này được gọi là bài toán đối ngẫu của bài toán ban đầu. Trong phần sau người ta 
sẽ chứng minh giá trị mục tiêu tối ưu của bài toán đối ngẫu bằng với giá trị mục tiêu tối 
ưu của bài toán gốc ban đầu. 

Người ta đưa bài toán đối ngẫu về dạng dể sử dụng bằng cách tính như sau : 

ø(y)= min { cIx+yT(b - Ax) })(x30) 

=min { cTx + yTb - yTAx }(x30) 

=min { yTb + (cT - yFA)x )(x30) 

= yTb + min { (cT - yTA)x } (x30) 

Ta thấy : 


min(c? — ? A)z = 
(z>0) 

0khi c” — A4 > 0 
Ụ 


[khôngxác đinh khi cˆ” — 7A < 0 


Vậy ta nhận được : 


ø(y) = yTb với cT-yTA 0 
Suy ra bài tóan đối ngẫu có dạng : 


maxg(w) = ”b 
UTA kế cT 
ụ€ R "tùy ý 


Hay là : 


maxg(0) = b”ụ 
ATu<c 
ục€ R "tùy ý 


Định nghĩa đối ngẫu trong trường hợp quy hoạch tổng quát 


Trong trường hợp quy hoạ ch tuyến tính tổng quát, những quy tắc sau đây được áp dụng 
để xây dựng bài toán đối ngẫu : 


- Hàm mục tiêu đối ngẫu : 

. max « min 

- Biến đối ngẫu : 

. Mỗi ràng buộc « một biến đối ngẫu 

- Chi phí đối ngẫu và giới hạn ràng buỘc : 
. Chi phí đối ngẫu « giới hạn ràng buộc 


- Ma trận ràng buộc đối ngẫu : 


. Ma trận chuyển vị 
- Chiều của ràng buộc và dấu của biến : 


. Ràng buộc trong bài toán max có dấu £ thì biến đối ngẫu trong bài toán min có dấu 3 0 ( 
trái chiều ) 


. Ràng buộc trong bài toán max có dấu = thì biến đối ngẫu trong bài toán min có dấu tùy ý. 


. Ràng buộc trong bài toán max có dấu 3 thì biến đối ngẫu trong bài toán min có dấu £ 0 ( 
trái chiều ) 


. Biến của bài toán max có dấu 3 0 thì ràng buộc đối ngẫu trong bài toán min có dấu 3 ( 
cùng chiều ) 


. Biến của bài toán max có dấu tùy ý thì ràng buộc đối ngẫu trong bài toán min có dấu = .. 


. Biến của bài toán max có dấu £ 0 thì ràng buộc trong bài toán đối ngẫu min có dấu £ ( 
cùng chiều ) 


Xét các ràng buộc dạng ma trận của một bài toán quy hoạch tuyến tính tổng quát như 
sau : 


đỚ11 G12... GỊi ... địn bị 
Ta 
T 
G; — Gị1 G2 Gịj Gin b¿ 
Số) 
(m1 đm2 ... Ômj ... Œmn In 
*„n 


Ký hiệu : 
œ7 là dòng thứ i(i=1,2,...,m) 
Aj là cột thứ j(=1,2,...,n) 


Khi đó, mối liên hệ giữa hai bài toán đối ngẫu có thể được trình bày như sau : 


Z(x) = cTx ® min 


xjJ30 
xjJ£0 


xj tự do 


Ví dụ 


w(y) = yTb ® max 


yi tự do 
vi£0 


yi30 


yTAj £ cj 
yTAj  cj 


yTAj = cj 


a- Hai bài toán sau đây là đối ngẫu : 


Ràng buộc / Dấu 


Cùng chiều 


Trái chiều 


max Z(%) = 30x + 10a 
2x +#a <4 
2XỊ + 2Xa < 6 


(P) 


%1;#2 báu 0 


{ 


mỉnw() = 4y + Öy› 
2y + 2y; > 30 
1+ 2ya > 10 


(D) 


U1;2 o 0 


{ 


b- Hai bài toán sau đây là đối ngẫu : 


minw(z) = #1 —Øa + a + 2X4 
#1 + 2xa — #4 + ðxa <€ 6 
2x1 — 3xa + 3xa — Á4xa > 7 
3XỊ = 2Xa - ĐX3 =ỹ 
fXxỊ †#as — 2xạ > 5 
(D) 


%1;#2 > 0,zztuy y ;#4 = 0 


tí 


maxz() = 6y + 7ya + 9ys + 5ya 
1ì + 2y¿ +3ys + 7y¿<1 
2ÿiT— 3ÿ¿T— 3ys ŠS —1 
Ý hi J J0 Poyg E4 
5y — 4y¿ — 2y¿ > 2 
(P) 


1i S 0,2 > 0,9atuy y,4 > 0 


t 


Đối với cặp bài toán đối ngẫu (P) và (D) chỉ xảy ra một trong ba trường hỢp sau : 
- Cả hai bài toán đều không có phương án tối ưu. 


- Cả hai bài toán đều có phương án, lúc đó chúng đều có phương án tối Ưu và giá trị hàm 
mục tiêu đối với hai phƯơng án tối ưu là bằng nhau. 


- Một trong hai bài toán không có phƯơng án, còn bài toán kia thì có phương án, khi đó bài 
toán có phương án không có phương án tối ưu. 


Các định lý về sự đối ngẫu 
a- Định lý 1 ( đối ngẫu yếu ) 


Xét hai bài toán đỐi ngẫu : 


(P) (Ð) 
max Z(#) = cTx min w(w) = bẨy 
Ấy ẤT >è 
z>0 y tùy ý 
tí tí 


Nếu # là phương án của bài toán (P) 
ÿ là phương án của bài toán (D) 


thì z() < (9) 


nghĩa là giá trị hàm mục tiêu của bài toán max không vượt quá giá trị hàm mục tiêu của bài 
toán đối ngầu min trên các phương án bất kỳ của mỗi bài toán . 


Chứng minh 


# là phương án của (P) nên : ÁZ = b 


by Az=p'b=b'y=u) 

g là phương án của (D) nên : 4ÿ > c 
bgTA>cT 

bợỤT AZ > cÝ# = z(®) 

Vậy z() < u(w) 


Định lý này được phát biểu và chứng minh cho hai bài toán đối ngẫu trong trường hợp 
tổng quát . 


b- Định lý 2 


Xét hai bài toán đỐi ngẫu : 


Œ) (0) 
max Z(#) = cTz min w() = by 
Ax=b ATfu>ec 
z>0 y tùy ý 


tí tí 
# là phương án khả thi của bài toán (P) 
ÿ là phương án khả thi của bài toán (D) 
Nếu z(#) = (2) thì z, y lần lượt là phương án tối ưu tương ứng của (P và (D). 
Chúng minh 
- Nếu # không là phương án tối ưu của bài toán (P) thì tồn tại một phương án x sao cho : 
Z(m) < z(#) 


b.(ÿ) < z(z) : điều này mâu thuẩn với định lý 1. 


- Nếu ÿ không là phương án tối ưu của bài toán (D) thì tồn tại một phƯơng án y sao cho : 
tu(w) < u() 

b.() < z(#) : điều này mâu thuẩn với định lý 1. 

Vậy z và lần lượt là phương án tối ưu của (P) và (D). 

c- Định lý 3 


Xét hai bài toán đối ngẫu : 


(P) (D) 
max Z(#) = cÏx min w() = bự 
Ash ATy>c 
z>0 y tùy ý 
tí tí 


Nếu x* là phương án tối ưu của bài toán (P) đối với cơ sở B thì phƯơng án tối ưu y* của 
bài toán (D) được tính bởi công thức : 


* 
Je th 
Chứng minh 


Do x* là phƯơng án tối ưu của (P) với cơ sở B nên thoả dấu hiệu tối ưu 
cẰ—c.B}A<0 
Bếp DA tực 


U* 


PTA >cT 


Py* là một phương án của (D) 
Mặt khác x* được tính bởi công thức : 
8u =3) DỊ (n0 


r¡gh 


và giá trị mục tiêu tối ưu của (P) là : 

Z(x*) = CÏX* = C1 

Ta có : 

00(UYabtu+ DU CSB l)' =(GbÐ n)0 
=0 n( 0) CbdipesEnitp Sở (ó | 


Theo định lý 2 thì y* là phương án tối Ưu của (D). 


Định lý này cho phép tìm phương án tối Ưu của bài toán quy hoạch tuyến tính đối ngẫu từ 
bài toán gốc. Trong đó : 


- c? được xác định trong bảng đơn hình tối ưu của (P). 


- B-1 gồm m cột tương ứng với m cột của ma trận cơ sở ban đầu lấy từ bảng đơn hình 
tối ưu của bài toán gốc. 


d- Định lý 4 ( sự đối ngẫu) 


Xét hai bài toán đối ngẫu 


(P) (D) 
max Z(#) = cTx min w(w) = bẨự 
Ax=b uc 
z>0 y tùy ý 
t tí 


- Nếu (P) và (D) đều có phƯơng án khả thi thì chúng có phương án tối ưu và giá trị của 
hàm mục tiêu tương ứng là bằng nhau. 


- Nếu một trong hai bài toán có phương án tối ưu không giới nội thì bài toán còn lại không 
có phƯơng án khả thi. 


Chứng minh 

- Đây là kết quả của định lý 3. 

- Giả sử rằng phương án tối ưu của (D) không giới nội, tức là tồn tại một phương án khả 
thi y của (D) sao cho w(y)= bTy nhỏ tuỳ ý. Điều này cũng có nghĩa là : với mọi M>0 lớn 
tuỳ ý luôn tìm được một phƯơng án khả thi của (D) sao cho : 

b ụ <—M 

Nếu (P) có phương án khả thi là # thì theo định lý 1 ta có : 

z(®) = c†% < u(ÿ) = b*Ụ < —M 

Điều này dẫn đến mâu thuẩn 

e- Định lý 5 (tính bổ sung ) 


Xét hai bài toán đối ngẫu 


Œ) (0) 
max Z(#) = cÏx min w() = bŸy 
AE ATgy>c 
z>0 y tùy ý 
tí tí 


#,„ là phương án khả thi tương ứng của (P) và (D). 


Điều kiện cần và đủ để Z,ÿ cũng là phƯơng án tối Ưu là : 


#7(ATg —- c1) =0 
Chứng minh 
- Do # là phương án khả thi của (P) nên : 
J1.—=U 
>S:( Am)? ii? 
= #TẠT — bĩ 
—> #T ATu — bĩụ 
=> #T AT — #øTc =— bTy-cTe(xfc = cTz) 
=> #T(ATg — c) = bTp-cT#(x) 


- Theo kết quả (*) : 


cxy — b Tụ 


. Nếu #,ÿ là phương án tối ưu của (P) và (D) thì theo định lý4_ => cfz — by =0 
=> #T7(ATg — c) =0 


. Nếu #?(ATg — e) = 0 > b#ụ — cfz = 0 => b*ụ = c†z 
Theo định lý 2 thì #,ÿ là phương án tối ưu. 
GIẢI THUẬT ĐỐI NGẪU 


Xét hai bài toán đỐi ngẫu : 


maxz(z) = cŸz minw(0) = bŸụ 
Ax=P APu Sẻ 
zœ>0 tuy 
(P) và (D) 
{ { 


Chúng ta sẽ xét xem giải thuật đơn hình cơ bản đã biết trong chương trước được áp dụng 
như thế nào đối với bài toán đối ngẫu. 


Giả sử rằng B là một cơ sở của bài toán (P) thoả : 


U=cCDH- và W ='ÉN 


#Ế#N — 0 
r¡gh 
Nếu B cũng là một cơ sở khả thi của bài toán gốc, tức là , thì (theo 
[ 
ng = 


định lý đối ngẫu) y, x lần lượt là phương án tối ưu của bài toán đối ngẫu và bài toán gốc. 
%B 
#N 
r¡gh 


Nếu không thì không là phương án của bài toán gốc vì zg = b = B”b không thể 


0. 
Để tiện việc trình bày ta xét (m=3 , n=5) : 


maxz(#) = C11 + Ca#2 -T C33 + CA%4 + CsZg 
Œ1121 + đ12#2 -E G133 + đ14#4 -E đ1s#s = ỦỊ 
đ211 + 0222 -E 233 + 0244 + gạs#s — Dạ 
G3121 -E đạ2#2 -E G333 -E G344 -E đạs2g —= Dạ 


(P) 


%1;12;T3;,14;#5 Boy 0 


tÍ 


Các dữ liệu của (P) đuợc trình bày trong bảng sau : 


x1 x2 x3 x4 x5 
c1 c2 c3 c4 c5 
a11 a12 a13 a14 a15 b1 
a21 a22 a23 a24 a25 b2 
a31 a32 a33 a34 a35 b3 


và bài toán đối ngẫu 


minw() = bi + ba› + Dạ13 
G1191 -F 6212 + 313 > CỊ 
G121 -F 0222 + 32s > Ca 
G131 -F 0232 + G333 > C3 
G141 -F 0242 + 0344 2 Cạ 
G1591 -F đ2592 T G353 > Cs 


(D) 


U1;Ù2 ;/stuy ở: 


tít 


Người ta đưa (D) về dạng chính tắc bằng cách thêm các biến phụ y4 y5, y6, y7,y8 0. 
Chúng không ảnh hưởng đến hàm mục tiêu. 


minw(0) = 191 + baa + bass + 0.0a + 0.0; + 0.0 + 0.0; + Ö.0g 
G111 -F G2192 -T G313 — 4 — C1 
G1291 -T G2292 -F G323 — Us — Ca 
G131 -F G2392 - G333 — 6 — C3 
0141 -T G2492 TT G344 — Uï — CẠ 
G151 -F đ2s12 - G353 — 1s — Cs 


LAI ;U2,stuy y -U4:U5;U6›:7;U5 bo 0 


tít 


Các dữ liệu của (D) được trình bày trong bảng sau : 


yl y2 y3 y4 y5 y6 y7 y8 

b1 b2 b3 0 0 0 0 0 

a11 a21 a31 sĩ 0 0 0 0 c1 
a12 a22 a32 0 sĩ 0 0 0 c2 
a13 a23 a33 0 0 S 0 0 c3 
a14 a24 a34 0 0 0 s1 0 c4 
a15 a25 a35 0 0 0 0 SỊ c5 


Giả sử rằng m cột đầu tiên của A là một cơ sở B của (P) thì hai bảng trên được trình bày 
rút gỌn nhƯ sau : 


B N b 


Bảng (P) 


vT y4....yÖ 
bT 0 
BT -Im 0 cB 


NT 0 -Ïn-m cN 


Bảng (D) 


Để đưa bài toán đối ngẫu về dạng chuẩn người ta nhân (bên trái) bảng (D) với bảng sau 
đầy : 


BNN -Ïn-m 


Khi đó người ta được bảng kết quả có dạng : 


ụ” _ y4y5y6 y7y8 
0 b=Bb 0 
m Im pg1T 0 2112 800 + 
Ân 0 CÔN V2 —B!N” In-m —CN = -c?. —c1B 1N” 


Bảng này cho ta một quy hoạch tuyến tính dạng chuẩn với ma trận đơn vị (cơ sở) tương 
ứng với các cỘt y1 y2 y3 y7 y8. 


Áp dụng giải thuật đơn hình cơ bản vào kết quả này cho ta quy tắc đổi cơ sở nhƯ sau : 
Tính :b = B †b> 0 
a- Nếu b > 0 thì giải thuật kết thúc, khi đó : 


ụ = cbÿB'' là phương án tối ưu của bài toán đối ngẫu . 


%B 
#N 


b 
0_ là phương án tối ưu của bài toán gốc . 


r¡gh 


b- Nếu tồn tại r sao cho „ € b,b„ < 0 thì xảy ra một trong hai trường hợp sau : 
- Nếu trong dòng r của Ñ có thành phần < 0 thì người ta tính : 


8 1 { li } 
—=— — IIIH4 — 
Ny Nụ 


VJ LÑ¡ <0 


Gi 


Như vậy : đối với bài toán đối ngẫu thì biến yr đi vào cơ sở và biến ys ra khỏi cơ sở, 
trong khi đó đối với bài toán gốc thì biến xs đi vào cơ sở và biến xr ra khỏi cơ sỞ. 
- Nếu mọi thành phần trong dòng r của Ñ đều > 0 thì phương án tối ưu của bài toán đối 


ngẫu là không giới nội, điều này (theo định lý đối ngẫu) dẫn đến bài toán gốc không có 
phương án. 


Ví dụ : Xét bài toán 


minw(#) = #1 — #3 
ø1— 2X T†#s = Ì 
#1 + 3xa + 4 = 2 


Bài toán đối ngẫu của (D) là : 


maxz(9) = 1ì + 2ya 
ú1i 9a S1 
—2y Tởya<0 
fM c 


2z<0 (®) 


tí 


y1, y2 là tùy ý 

Ta có thể chọn bài toán (D) hoặc (P) để giải tìm phương án tối ưu bằng phương pháp đơn 
hình, từ đó suy ra phƯơng án tối ưu của bài toán còn lại theo kết quả trên. Trong ví dụ này 
ta chọn bài toán (D) để giải vì có chứa sẵn ma trận đƠn vị. 


Giải bài toán (D) bằng phương pháp đơn hình cải tiến ta được : 


sỉ 2 -2 0 
CB, #Py #1 %2 %3 
sỉ 3 Ÿ 0 1 
0 2 3 T 0 
cT Ỷ 0 -1 
sử " 0 0 


Giải thuật dừng vì thoả dấu hiệu tối ưu của bài toán min. 


Phương án tối ưu của bài toán (D) là : 


Suy ra phương án tối ưu của (P) là : 


1 2+ 
u? =l0i 9a| =cpB~ a0 ¡ =lL Số 
| 
sb)=##y=llÍ ".ủ 
{ 


CÂU HỎI CHƯƠNG 3 


Z4 


©a|k2 


w(%0) 


1- Bạn hiểu như thế nào về khái niệm đối ngẫu ? 


2- Quy hoạch tuyến tính đối ngẫu của một quy hoach tuyến tính chính tắc có dạng như 
thế nào ? 


3- Bạn hãy nêu ra các quy tắc đối ngẫu. Cho ví dụ. 


4- Giá trị hàm mục tiêu của hai quy hoạch tuyến tính đối ngẫu thì như thế nào ? . Chứng 
minh 


BÀI TẬP CHƯƠNG 3 

1- Xét bài toán quy hoạch tuyến tính 

max Z = 7x1 + 5x2 

2X, OXC. lƠ 

(P)2x1+x2 13 

2X. to 

3x1 18 

x1... x2“ .0 

a- Tìm bài toán đối ngẫu (D) từ bài toán (P) 
b- Tìm phương án tối ưu cho bài toán (P) 

c- Từ bảng đơn hình tối ưu của (P). Hãy tìm phương án tối ưu cho bài toán (D) 
2- Xét bài toán quy hoạch tuyến tính 

min w= x1 + x2 

x1 -2x3+x4=2 

(D)x2 - x3 + 2x4 = 1 

x3-x4+x5=5 

xI 0, 1=15 

a- Tìm bài toán đối ngẫu của bài toán (D) 


b- Tìm phương án tối ưu của bài toán (D) 


c- TỪ bảng đơn hình tối ưu của bài toán (D). Hãy tìm phƯơng án tối ưu cho bài toán đối 
ngẫu Ở câu a. 


3- Xét bài toán quy hoạch tuyến tính 
min w = -2x1 - x4 

xi +x2+5x535=20 

(D)x2+2x4 5 

XI +x2sx3 8 

xi tùy ý (¡=1 4) 


Tìm bài toán đối ngẫu (P) của bài toán (D). Từ bài toán (P) hãy chỉ ra rằng (P) không tồn 
tại phương án tối ưu do đó (D) cũng tồn tại phương án tối ưu. 


4- Cho bài toán quy hoạch tuyến tính 


maxz = 2xị + 4xaạ + #3 + Ø4 
#+ +3xạ +z#a <€ 1 
—ÖXa — 2xa < 3 
4xa + 4xs + xạ < 3 


(D) #ï.:0(7=1s 4) 


1- Tìm bài toán đối ngẫu của bài toán đã cho. 
2- Giải bài toán đã cho rồi suy ra kết quả của bài toán đối ngẫu. 


5- Cho bài toán quy hoạch tuyến tính 


maxz = 27x + 50Za + l8zs 
#1 + 2xạ +za <S 2 
—2xI + #ạ — 2xạ < 4 
#1 + 2xạ — 4x < —2 
#1,zatuú ý, xạ € Ủỳ 


(D) 


t 


a- Tìm bài toán đối ngẫu của bài toán đã cho. 


b- Giải bài toán đối ngẫu rồi suy ra kết quả của bài toán đã cho. 


Ứng dụng quy hoạch tuyến tính 

Chương này trình bày các bài toán để thấy khả năng ứng dụng rộng rãi của quy hoạch 
tuyến tính. Bài toán trò chơi được trình bày một cách chỉ tiết, các bày toán còn lại chỉ trình 
bày mô hình. Việc giải các bài toán này được nghiên cứu thêm trong các môn tiếp theo. 


MỞ ĐẦU 

Trong thực tế hay gặp tình huống là phải chọn một quyết định (bấp bênh) do phải đối mặt 
với một đối thủ thông minh và có quyền lợi đối lập với ta : ví dỤ trong các trò chơi tranh 
chấp, trong quân sự, trong vận động tranh cử... 

Nghiên cứu việc chọn quyết định trong những trường hợp đối kháng này có tên gọi là lý 
thuyết trò chơi. Ở đây người chọn quyết định và đối thủ đều được gọi là người chơi. Mỗi 
người chơi có một tập hợp các hành động để lựa chọn được gọi là chiến lược. 


Chúng ta xét mỘt trường hợp đơn giản là trò chơi hai người : phần thưởng sẽ là cái được 
của một người và chính là cái mất của người kia. 


Giải một trò chơi nghĩa là tìm chiến lược tốt nhất cho mỗi người chơi. Hai người chơi 
thường được ký hiệu là A và B, chiến lược tương ứng của mỗi người được ký hiệu là : 


A:i(=1m) 
B:j(=1n) 


Giải thưởng ứng với chiến lược (¡,j) của hai người được ký hiệu là aij và được viết thành 
một bảng như sau : 


b 1 2 n 

A 

1 a11 a12 a1n 
2 a21 a22 vã a2n 
m am1 am2 Sà amn 


Ví dụ : 


Ầ 1 1 0 9) 1 
9) 2 2 1 0 
ki -1 xỉ 0 3 
Đối với A : 


- Nếu A đi nước 1 (dòng 1) thì A sẽ : 

. Thắng 1 điểm nếu B đi nước 1(thắng) 

. Thắng 0 điểm nếu B đi nước 2(hoà) 

. Thắng -2 điểm nếu B đi nước 3(thua) 

. Thắng 1 điểm nếu B đi nước 4(thắng) 

Những trường hợp còn lại là tương tự . 

Đối với B : 

- Nếu B đi nước 2 (cột 2) thì B sẽ : 

. Thua 0 điểm nếu A đi nước 1 

. Thua 2 điểm nếu A đi nước 2 

. Thua -1 điểm nếu A đi nước 3 

Những trường hợp còn lại là tương tự. 

Ñ ghiệm tối ưu của trò chơi, có khi gọi tắt là nghiệm, là bộ chiến lược (¡*,j*) có tính chất là 
nếu một người lấy chiến lược khác còn người kia vẫn giữ nguyên thì phần thưởng cho 


người đi khác sẽ bị thiệt hại. Giải trò chơi có nghĩa là tìm nghiệm tối Ưu. 


BÀI TOÁN TRÒ CHƠI 


Trò chơi có nghiệm ốn định 

Hai nhà chính trị A và B vận động tranh cử 1 ghế ở nghị viện trong 2 ngày cuối quan trọng 
nhất ở hai thành phố P và Q. Mỗi người phải đặt kế hoạch vận động mà không biết được 
kế hoạch của đối phương. Các cố vấn đưa ra 3 chiến lược : 

- Ở mỗi thành phố một ngày 


- Ở cả 2 ngày ở thành phố P 


- Ở cả 2 ngày ở thành phố Q và đánh giá kết quả vận động tương Ứng như sau : 


1 2 3 B 
A 1 1 b) 4 

2 1 0 5 

3 0 1 v† 


Dữ liệu là tổng số phiếu, tính theo đơn vị là ngàn, mà A sẽ dành được từ B hay ngược lại . 


Đây là một trường hợp đơn giản mà người ta có thể giải được bằng khái niệm chiến lược 
bị trội hơn như sau : 


- Đối với A thì chiến lược 3 bị trội hơn bởi chiến lược 1 và 2 vì nó mang đến cho A số 
điểm thắng ít, nên dù B có chọn chiến lược nào thì A cũng vẫn chọn chiến luợc 1 hoặc 2 
mà bỏ chiến lược 3. Ta có : 


A 1 ƒ 2 Á 
2 1 0 5 
3 0 1 -† 


- Đối với B thì chiến lược 3 bị trội hơn bởi chiến lược 1 và 2 vì nó mang đến cho B số 
điểm thua nhiều nên B bỏ chiến lược 3. Ta có : 


1 Ũ 3 B 
A 1 1 2 Á 

2 ƒ 0 5 

3 0 T Ỉ 


- Đối với A thì chiến lược 2 bị trội hơn bởi chiến lược 1 vì vậy A bỏ chiến lược 2. Ta có : 


A 1 ƒ 2 Á 
2 1 0 5 
3 0 1 s† 


- Đối với B thì chiến lược 2 bị trội hơn bởi chiến lược 1 vì vậy B bỏ chiến lược 2. Ta có : 


1 2 3 B 
A 1 1 ? Á 
2 1 0 5 
3 0 1 +1 


Cuối cùng thì bộ chiến lược (1,1) là nghiệm tối ưu của trò chơi với kết quả là người A thu 
thêm được 1 (ngàn) phiếu từ người B. 


Trong nhiều trường hợp, khi dùng chiến lược bị trội hơn chỉ mới giảm được cở của bài 
toán mà chưa giải quyết xong vấn đề đặt ra. 


Chiến lược MaxiMin và MiniMax 


Xét ví dụ tương tự như ví dụ trên nhưng bảng kết quả vận động được các cố vấn đánh giá 
nhƯ sau : 


1 5; Đ) B 
A 1 -3 K) 6 
Bộ 1 0 2 
3 5 s2 -4 


Đây là trường hợp người chọn quyết định nghĩ là đối phương thông minh và cố ý chọn 
quyết định chống lại mình nên họ luôn nghĩ đến chiến lượt “ăn chắc”, đó là MaxiMin(A) và 
MiniMax(B) nhƯ sau : 

a- MaxiMin(A) 


A luôn xem B là đối thủ thông minh. Khi A đi nước ¡0 (dòng ¡0) thì B sẽ chọn nước đi j0 
(cột j0) sao cho A thắng điểm ít nhất . Nghĩa là B đi vào ô : 


đjpjạ Mintaij} 


Trong tình huống đó A sẽ chọn nước đi sao cho A thắng nhiều điểm nhất. Chiến thuật của 
A là đi vào ô : 


L 7 
A đi nước 1 thì B sẽ đi nước 1: a11=-3 


A đi nước 2 thì B sẽ đi nước 2 : a22=0 


A đi nước 3 thì B sẽ đi nước 3 : a33=-4 


A 1 -3 m9) 6 
5 1 0 Mu 
3 5 -Z -4 


Vậy MaxiMin(A) = a22 = 0 
b- MiniMax(B) 


B luôn xem A là đối thủ thông minh. Khi B đi nước j0 (cột j0) thì A sẽ chọn nƯớc đi ¡0 
(dòng i0) sao cho B thua điểm nhiều nhất . Nghĩa là A đi vào ô 


Ôjsjy — max{aj, } 
V¿ 


Trong tình huống đó B sẽ chọn nước đi sao cho B thua ít điểm nhất. Chiến thuật của B là đi 
Vào Ô : 


9p — ajy — MiniMax(P) = mìn | max{aạ}Ì 


B đi nước 1 thì A sẽ đi nước 3: a31=5 

B đi nước 2 thì A sẽ đi nước 2 : a22=0 

B đi nước 3 thì B sẽ đi nước 1 : a13=6 

Vậy MiniMax(B) = a22= 0 

Lần này ta thấy rằng : 

MaxiMin(A) = MiniMax(B) = a22= 0 

Bộ chiến lược (2,2) có giá trị là 0 là nghiệm tối ưu của trò chơi vì nếu người nào đi lệch và 


người kia đi đúng thì người đi đúng thu lợi nhiều hơn giá trị của trò chơi. Nghiệm tối ưu 
trong trường hợp này còn được gọi là nghiệm ổn định. 


Trò chơi không có nghiệm không Ổn định 


Xét ví dụ tương tự như trên với bảng kết quả được các chuyên gia đánh giá nhƯ sau : 


Khi A và B dùng chiến lược MaxiMin và MiniMax của mình thì cho kết quả như sau : 
MaxiMin(A) = a12 = -2 
MiniMax(B) = a13 = 2 


Vì MaxiMin(A) và MiniMax(P) là khác nhau nên trò chơi không có nghiệm ổn định. Ta xem 
điều gì có thể xảy ra ? 


- A tính rằng nếu B thực hiện đúng chiến lược của mình là chọn cột 3 thì A sẽ chọn chiến 
lược 1 để thắng 2 từ B (thay vì thắng -2) 


1 2 3 B 
A ï 0 “7 2 

5) 5 4 sử 

3 2 3 -4 


- Lúc này B sẽ suy tính và thấy rằng phải chọn chiến lược 2 để thua -2 từ A (thay vì thua 2). 


- Đến lượt A cũng đủ thông minh để tính liền được 2 nước, biết được B sẽ chọn chiến 
lược 2 nên A sẽ dùng chiến lược 2 để thắng 4 từ B. 


1 2 3 B 
A † 0 d2 2 

2 5 4 -3 

Ầ 2 3 -4 


- Nhưng B cũng tính được điều này nên sẽ quay lại chọn chiến lược 3 để thua -3 từ A.. 


- Cũng như B, A cũng sẽ tính được điều này nên sẽ quay lại chọn chiến lược 1 để thắng 2 
từ B. 


1 2 3 B 
Ầ l| 0 su, 2 

5 5 4 ^ủ 

3 ? 3 -4 


Như vậy ta đã xoay đúng một vòng, và nếu cứ lập luận như vậy thì ta sễ xoay vòng mãi. 
Những bộ chiến lược nhận được trong khi xoay vòng là những nghiệm không ổ định. 


Chiến lược hỗn hợp 

Để có được lời giải của trò chơi không có nghiệm ổn định người ta đưa ra khái niệm chiến 
lược hỗn hợp. Mỗi người chơi không chọn một chiến lược thuần túy như trước đây mà 
chọn một phân bố xác suất sử dụng tất cả các chiến lược. 


Xét trò chơi giữa A và B có ma trận điểm dƯơng có dạng tổng quát : 


ái Ớ1]1 G12 
A 

2 G21 G22 

In G1 đn 2 


Giả sử rằng : 

MaxiMin(4) = a¿„;¿ = gA 

MiniMax(Ð) = a;¡„;, = gp 

G2xj4 Z” đáyjp 

GỌI : 

.pi>0(=1 m) là tần suất nước đi thứ ¡ của A với 
p1 +p2+... + pm = 1 

.qj>0(=1 n) là tẫn suất nước đi thứ j của B với 


q1 + q2+...+ qn= 1 


q1 q2 


ỚIn 


đ2n 


Cmn 


qn 


p1 1 Ớ11 G12 si đIn 


p2 2 G21 G22 Net G2n 


pm m đựy¡ï đựn2 vác Địa 


Vấn đề đặt ra là : 


-Tìm tần suất pi > 0 của nước đi thứ ¡ (¡ =1 m) của A sao cho đối với mỗi nước đi thứ j 
của B số điểm thắng trung bình của A không nhỏ thua gA : 


p1a1j + p2a2J +..... + pmamJ( J= 1 n) 

Cũng có nghĩa là tìm pi sao cho : 

p1a1j + p2a2J +..... +pmamj g1 gA(j=1 n) 
g1 max 


- Tìm tần suất qj > 0 của nước đi thứ j (j =1 n) của B sao cho đối với mỗi nước đi thứ ¡ của 
A số điểm thua trung bình của B không lớn hơn gB : 


qlai1 + q2al2 +.... + qnain(Ii=1 m) 

Cũng có nghĩa là tìm các qj sao cho : 

q1ail + q2ai2 + ..... +tanan gø2 gB(1=1 m) 
ø2 min 


Khi đó hai bài toán quy hoạch tuyến tính thu được là : 


InaXg1 min 
01+ Ø2 +... + Ðm = Ì 
i0 0388 EssBmdag S1 =>8) 


ø¿ > 0(0 = 1 — n) 


t 


minga max 


01 +©03+‹.tu, <1 
g; > 0 =1 ->m) 


t 


Chia các ràng buộc của bài toán thứ nhất cho g1>0 và đặt : 
#¡ = Ÿ. ( = 1 — Tn) 

Chia các ràng buộc của bài toán thứ hai cho ø2>0 và đặt : 
U¿ = „(j= 1— n) 

Khi đó hai bài toán quy hoạch tuyến tính trên trở thành : 


mỉnz- =#1 +#as +... + mm 
đ1j#1 + đ2j2 +... + đmj#m > 1(7 =l1— m) 
ø; > 0(2—= 1 >mn 
(Ð) ( ) 


t 


max =1 T12 +... 13 
;>0(7—=1—n 
(P) Ùj ( ) 


tí 


Đây là hai bài toán đối ngẫu . Chọn một trong hai để giải 
Ví dụ : 


Xét trò chơi giữa A và B có bảng điểm như sau : 


A 1 -1 2 
2 1 z2 
3 3 4 


Theo chiến thuật của A và của B ta có : 
MaxiMin(A) = a11 
MiniMax(B) = a23 


Tăng đồng loạt các ô của bảng điểm lên 4 ta được : 


Gọi 


pi 0 là tần suất nước đi thứ ¡ của A (¡=1 3) 


pl+p2+p3=1 


qj 0 là tần suất nước đi thứ j của B (=1 3) 


q1 + q2 + q3 =1 


Thực hiện tương tự như trên ta được hai bài toán đối ngẫu như sau : 


q1 


p2 


p3 


minWw = j- =đi +#¿ +#s 
JxI + ðxa + 7xa > Ì 
6xi + 2xa + 8xa > 1 
ð5xị +xa +#s> 
Xr> 0;x¿ >0xa >0 


(D) 


t 


Ta chọn bài toán (P) để giải. 


Đưa bài toán (P) về dạng chuẩn : 


(P) 


q2 q3 
&) 6 5 
5 2 6 
7 8 Ï 
maxz = - —= 1 +12 t3 


0 
3y¡ +ổy¿ +ðy; SÌ 
Syi †2ya †Ôya <1 
fyị t8ya + 1s S1 
ú1i>0,yạ >0,yạ>0 


t 


max øZ = : = 1i +12 + 0s + 0:4 + 0-0 + Ö-U& 


(P) 


đyi 
3y] 


Õy; 


5Ys +14 = 


2y; 


6y¿ + s = Ì 


fyị tổỔy¿+ 13 +T 1s = ] 
j >0;ys>0;ys >0,Yy¿>0,ys S50;ye>0 


tiit 


Dùng giải thuật đơn hình cải tiến : 


Cbạ 


CB;ị 


?Bụ 


?Bì 


LAI 


Lai 


12 L) U4 Uö 1ö 

6 5 1 0 0 

) 6 0 1 0 

8 ll 0 0 1 

1 1 0 0 0 

lÌ 1 0 0 0 
12 L) U4 U5 1ö 
# # |: |°o |-‡ 
-# |# |° |: |- 
P1 0 09 
1 1 0 0 0 


Z0 


| 


Z1 


sĩ ni, 
CB; 1B; LI 2 

0 4 0 “S 
H 3 0 —Ÿ* 
1 1 1 mn 
cP 1 1 

Z 0 3 
CB; 1B; 1 2 

ƒ 2 0 1 

1 3 0 0 

1 1 1 0 

cT 1 1 

G. 0 0 


Phương án tối Ưu của bài toán (P) là : 


U3 


KE) 


¬|= 


U4 


— 214 


| 


| 


SUY ra 


Phương án tối ưu của bài toán đối ngẫu (D) được tính bằng công thức sau : 
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j7 16 7 

214 107 214 

LING — ,Tn-1 _—_ 13 9 12 
m° |@t,.đa i8ú| S=£@nÐ “ez|1.:] J| ca ƒnỸớ. # in: 
“. .. .ằẮ. 

107 T07 107 


=l# ##. #r] 


214 107 214 


sí 
214 
HÙ T 10 ..-: 


K 
214 


t 


BÀI TOÁN VẬN TẢI 


Mở đầu 


Bài toán vận tải là bài toán quan trọng nhất trong các bài toán quy hoạch tuyến tính. Người 
ta tổng kết rằng 85% các bài toán quy hoạch tuyến tính gặp trong ứng dụng là bài toán vận 
tải hoặc mở rộng của nó. Thuật ngữ bài toán vận tải thường được hiểu là bài toán vận 
chuyển sao cho cước phí nhỏ nhất. 


Các khái niệm cơ bản 


Bài toán vận tải được mô tả như là một bài toán về dòng dữ liệu gồm tập hợp các nút N 
được chia thành hai phần rời nhau : các nút nguồn S và các nút đích D, tức là : 


WN=S5UÐĐ 
snD=Jỹ 


và mỗi cung (i,j) trong tập các cung A đều có gốc trong S và có ngỌn trong D. 


S:Các nút nguồnD:Các nút đích 


Các nút thuộc S được gọi là các nút nguồn (cung), các nút thuộc D được gọi là các nút đích 
(cầu). Một cách tổng quát, bài toán vận tải trình bày được bằng đồ thị. 


Ở bài toán vận tải đôi khi còn có thêm giả thiết nữa là mỗi nút nguồn đều có cung nối với 
mọi nút đích. Ở đây ta chỉ đề cập đến bài toán vận tải có thêm giả thiết này và sẽ gọi tắt là 
bài toán vận tải. 

Đối với bài toán vận tải người ta thường ký hiệu 

sỉ S là nguồn phát Ở nút i(=1 m) 

dj_ D là nhu cầu thu của nút j (=1 n) 

Trong trường hợp các nguồn phát không chuyển hết sang các nút cầu vì đã đủ nhu cầu thì 
bài toán vận tải được gọi là bài toán vận tải mở. Có thể đưa một bài toán vận tải mở về 
một bài toán vận tải (đóng) bằng cách thêm vào một nút cầu giả thứ (n+1) với nhu cầu 
được xác định nhƯ sau : 


đu+i = 3 34< 8i — )Jj— đỳ 

Bài toán vận tải cân bằng thu phát 

a- Thiết lập bài toán 

Có m nơi A1, A2,...,Am cung cấp một loại hàng với khối lượng tương ứng là a1, a2,....,am. 
Hàng được cung cấp cho n nơi B1, B2,...., Bn với khối lượng tiêu thụ tương ứng là b1, 
b2,....,bn. 

Cước phí chuyên chở một đơn vị hàng từ điểm phát Ai đến điểm thu Bị là cij. 

Hãy lập kế hoạch vận chuyển từ mỗi điểm phát đến mỗi điểm thu bao nhiêu hàng để : 
- Các điểm phát đều phát hết hàng 

- Các điểm thu đều nhận đủ hàng 

- Tổng cước phí phải trả là ít nhất 

Gọi xij là lượng hàng chuyển từ điểm phát Ai đến điểm thu Bj,xij 0. 


Vì tổng lượng hàng phát đi từ mỗi điểm phát Ai đến mọi điểm thu Bj bằng lượng hàng phát 
tỪ Ai nên : 


#¿i - ¿2 +.... T† #ịn —= d;(¿ == 1;2;...;0) 


Vì tổng lượng hàng thu được tại mỗi điểm thu Bịj từ mọi điểm phát Ai bằng lượng hàng 
cần thu tại Bj nên : 


#1j + 223j +.... T #mj — b;¡(J D0 n) 
Để tổng cước phí là ít nhất cần phải có : 


min2z(#) = ồ}; By Cij3ij 


Với các phân tích trên ta có mô hình của bài toán nhƯ sau : 


mina(ø) = 37 327 ¡ayey() 
Suy g + — Hi(U =2 a21n) 
3, 18y=07=12. 10) 
(2) 


Phương án - PhƯơng án tối ưu 


Một ma trận X=[xij]m.n thỏa (2) và (3) được gọi là phương án, thỏa thêm (1) được gọi là 


phương án tối ưu. 


b- Dạng bảng của bài toán vận tải 


Có thể giải bài toán vận tải theo cách của quy hoạch tuyến tính. Tuy nhiên do tính chất đặc 
biệt của bài toán vận tải nên người ta nghĩ ra một thuật toán hiệu quả hơn. Trước tiên 
người ta trình bày bài toán vận tải dưới dạng bảng như sau : 


ThuCướcPhát b1 b2 

a1 c11x11 c12x12 ........ 
a2 c21x21 È22x22ˆ |‹‹...‹: 
ai cI1xI1 c2x12  |........ 


bj 
c1jx1j 


c2jx2j 


cijxij 


bn 


clnx1n 


c2nx2n 


cinxin 


am cm1xm1 cm2xm2 | ........ cmjxmj | ........ cmnxmn 
Trong bảng mỗi hàng mô tả một điểm phát, mỗi cột mô tả một điểm thu, mỗi ô mô tả một 
tuyến đường đi từ một điểm phát tới một điểm thu. 
Dây chuyền - Chu trình 
Một dãy các ô của bảng mà hai ô liên tiếp nằm trong cùng một hàng hoặc một cột, ba ô liên 


tiếp không cùng nằm trên một hàng hoặc một cột được gọi là một dây chuyền. Ta thấy 
rằng hai ô liền nhau trong một dây chuyền có chỉ số hàng hoặc chỉ số cột bằng nhau 


Dây chuyền : (1,2) (1,3) (2,3) (2,4) (4.4) (4,1) 


Một dây chuyền khép kín, ô đầu tiên và ô cuối cùng bằng nhau, được gọi là một chu trình.Ta 
thấy rằng số ô trong một chu trình là một số chẵn. 


Chu trình : (1,1) (1,3) (2,3) (2,4) (4,4) (4,1) (1,1) 


Ô chọn - Ô loại 


Giả sử ma trận X=[xij]m.n (i=1,2,...,m) (j=1,2,...,n) là một phương án của bài toán vận tải. 


Những ô trong bảng tương ứỨng với xij >0 được gọi là ô chọn, những ô còn lại được gọi là ô 
loại. 


Phương án cơ bản 


Một phương án mà các ô chọn không tạo thành một chu trình được gọi là phương án cƠ 
bản. 


Một phương án có đủ m+n-1 ô chọn được gọi là không suy biẾn, có ít hơn m+n-1 ô chọn 
được gọi là suy biến. Trong trường hợp suy biến người ta chọn bổ sung vào phương án cơ 
bản một số ô loại có lượng hàng bằng 0 để phương án cơ bản trở thành không suy biến 

c- Giải bài toán vận tải 


Xét bài toán vận tải có số lượng phát, số lượng thu và ma trân cước phí ở dạng bảng như 
Sau : 


30 20 60 
90 Đ 4 1 
40 3 2 6 
70 ⁄ 9 IT 


LẬP PHƯƠNG ÁN CƠ BẢN BAN ĐẦU 


Phương án cơ bản ban đầu được xác định bằng cách ưu tiên phân phối nhiều nhất vào ô có 
cước phí nhỏ nhất (r,s) ( gọi là ô chọn). Khi đó : nếu điểm phát r đã phát hết hàng thì xóa 
hàng r của bảng và số lượng cần thu tại điểm s chỉ còn là bs-ar ; nếu điểm thu s đã nhận đủ 
hàng thì xóa cỘt s của bảng và số lượng phát còn lại tại điểm phát r là ar-bs 


Bảng mới thu được có kích thước giảm đi. Tiếp tục phân phối như trên cho đến khi hết 
hàng. 


Các ô chọn trong quá trình phân phối, sẽ không chứa chu trình, là một phương án cơ bản. 


Nếu phương án cơ bản suy biến, chưa đủ m+n-1 ô, thì bổ sung thêm một số " ô chọn 0 " 


Áp dụng vào bài toán đang xét : 


1- Phân vào ô (1,3) 50 .. Hàng (1) bị xóa . Cột (3) còn thu 60-50=10 


g0 20 
0 Đ 4 
40 3 2 
70 + 9 


2- Phân vào ô (2,2) 20. Cột (2) bị xóa . Hàng (2) còn phát 40-20=20 


30 0 
0 5 4 
20 3 220 
70 G 9 


3- Phân vào ô (2,1) 20.. Hàng (2) bị xóa . Cột (1) còn thu 80-20=60 


60 0 
0 5 4 
0 320 220 
70 7 9 


4- Phân vào ô (3,1) 60 . Cột (1) bị xóa . Hàng (3) còn phát 70-60=10 


0 5 4 150 
0 320 220 6 
10 760 9 11 


5- Phân vào ô (3,3) 10. Hết hàng. 


0 0 0 
0 5 4 150 
0 320 220 6 
0 760 9 1110 


Đã có 5 ô được chọn, chúng tạo thành một phương án cơ bản không suy biến vì số ô bằng 
với m+n-1=3+3-1. 


THUẬT TOÁN "QUY 0 CƯỚC PHÍ CÁC Ô CHỌN" 

Định lý 

Nếu cộng vào hàng ¡ và cột j của ma trận cước phí C=[cij] một số tùy ý ri và sj thì bài toán 
vận tải mới với ma trận cước phí mới C'=[c'ij=cij+ri+sj] thì phương án tối ưu của bài toán 
này cũng là phương án tối ưu của bài toán kia và ngược lại. 


Thuật toán "Quy 0 cước phí các ô chọn" gồm ba giai đoạn. 


Giai đoạn 1 : Quy 0 cước phí các ô chọn 


Sau khi xác định được phương án cơ bản có m+n-1 ô chọn, người ta cộng vào mỗi hàng ¡ và 
mỗi cột j của ma trận cước phí C=[cij] một sỐ ri và sj sao cho ma trận cước phí mới C' tại 
các ô chọn thỏa c1J=cij+ri+sj=0. 


Tiếp tục ví dụ trên ta thấy : 


320 220 6 r2=0 
760 9 11 10 r3=-4 
s1=-3 s2=-2 s3=-7 


Các giá trị cỘng vào phải thỏa hệ phương trình : 


l+ri+sa=0 
3+ra+sị=0 
2-+rsa+saạ=0 
f+raza+sị=0 


1l+rạ+sa=0 


thiệt 


Chọn r2=0, giải hệ ta được kết quả trên 


Ma trận cước phí mới thu được là : 


ö ö 0 50 
020 020 +1 
0 60 3 010 


Giai đoạn 2 : Kiểm tra tính tối ưu 


Sau khi quy 0 cước phí các ô chọn nết : các ô loại đều có cước phí 0 thì phương án đang 
xét là tỐi ưu, ngược lại thì chuyển sang giai đoạn 3 


Trong ví dụ này ta chuyển sang giai đoạn 3. 


Giai đoạn 3 : Xây dựng phương án mới tốt hơn 

1- Tìm ô đưa vào. 

Ô đưa vào là ô loại (¡*,j*) có cước phí nhỏ nhất và trở thành ô chọn 
Trong ví dụ này là ô (2,3). 

2- Tìm chu trình điều chỉnh. 


Chu trình điều chỉnh được tìm bằng cách bổ sung ô (¡*,j*) vào m+n-1 ô chọn ban đầu, khi đó 
sẽ xuất hiện một chu trình duy nhất, gọi là chu trình điều chỉnh V . 


Trong ví dụ này chu trình điều chỉnh là : 
Ÿ;0/3)(35)3.1)(2/1)025) 
3- Phân ô chẵn lẻ cho chu trình điều chỉnh. 


Đánh số thứ tự các ô trong chu trình điều chỉnh V bắt đầu từ ô (¡*,j*). Khi đó chu trình điều 
chỉnh V được phân thành hai lớp : 


VC : các ô có số thứ tự chẵn. 
VL : các ô có số thứ tự lẻ. 
4- Tìm ô đưa ra và lượng điều chỉnh. 


Trong số các ô có thứ tự chẵn chọn ô (r,s) được phân phối ít hàng nhất làm ô đưa ra, trở 
thành ô loại. Lượng hàng xrs Ởở ô đưa ra gọi là lượng điều chỉnh. 


Trong ví dụ này ô đưa ra là ô (3,3), lượng điều chỉnh là 10. 
5- Lập phương án mới. 


Phương án mới có được bằng cách thêm hoặc bớt lượng điều chỉnh trên chu trình điều 
chỉnh nhƯ sau : 


Ô có thứ tự chẵn bị bớt đi lượng điều chỉnh. 

Ô có thứ tự lẻ được cộng thêm lượng điều chỉnh. 

Ô ngoài chu trình điều chỉnh không thay đổi 

Trong ví dụ này ta thấy những ô trong chu trình điều chỉnh có sự thay đổi nhƯ sau : 
Ô (2,3) được thêm 10 trở thành 10 


Ô (3,3) bị bớt 10 trở thành 0 


Ô (3,1) được thêm 10 trở thành 70 


Ô (2,1) bị bớt 10 nên trở thành 10 


Khi đó phương án mới là : 
lô lô 
010 020 
070 3 


Quay về giai đoạn 1. 


Giai đoạn 1 : Quy 0 cước phí ô chọn 


lô lô 0 50 

010 020 -1 10 

070 3 0 

s1=0 s2=0 s3=1 
Ma trận cước phí mới là : 

7 „, 


010 0 20 


050 


010 


070 3 


Giai đoạn 2 : Kiểm tra tính tối ưu 


Đây là phƯơng án tối ưu 


30 20 
90 5 4 
40 3 10 220 
70 770 9 


Với cước phí là : 
1.50+3.10+2.20+6.10+7.70=670 


Khi sử dụng phương án ban đầu 


30 20 
90 Đ 4 
40 320 2.0 
70 760 9 


thì cước phí là : 
1.50+3.20+2.20+7.60+11.10=680 


Các bài toán được đưa về bài toán vận tải 


60 


150 


610 


11 


Có nhiều bài toán thực tế có tính chất không phải là ”vận tải ” nhưng có mô hình toán học 
là bài toán vận tải. Một số bài toán như vậy là : 


a- Bài toán bổ nhiệm 
Giả sử tập hợp S gồm m người và tập hợp D gồm n công việc (chức vụ). Cước phí của 
việc bổ nhiệm người ¡ S vào việc j D là cij (=1 m, j=1 n). Bài toán đặt ra là tìm cách chia 
mỗi người đúng một việc sao cho cước phí bổ nhiệm là nhỏ nhất. 
Người ta đặt biến (biến trên dòng) như sau : 
1nÕu ngêi ¡ nhẼn viÖc j 

0nÕu trêng hp kh c 


®%j = { 


thì bài toán trở thành : 

min È ”;-s È ,;cp Cij2ij 

Vì mỗi người nhận đúng 1 việc nên : 3 ;ep # = 1(W¿ € 6) 

Vì mỗi việc chỈ giao cho mỘt người nên : 3”,_s z¡„ = 1(V7 € D) 

Đây là bài toán vận tải nhưng có thêm yêu cầu là các biến xij chỉ lấy giá trị 0 hoặc 1. 

Bài toán bổ nhiệm cũng có khi được gọi là bài toán chọn (Choice Problem). Nhiều bài toán 
thực tế đa dạng có mô hình toán học là bài toán bổ nhiệm, chẳng hạn như bài toán phân bố 
hoả lực vào mục tiêu cần tiêu diệt. 

b- Bài toán vận tải với cung ít hơn cầu 


Xét một bài toán một bài toán vận tải với S là tập hợp m nút cung và D là tập hợp n nút cầu 
mà tổng nguồn cung nhỏ hơn tổng nhu cầu, tức là 


be, 8< ).a đ; 


Trong trường hợp này tất nhiên không thể đáp ứng đủ nhu cầu dj cho mỗi nút j=1 n cho nên 
ràng buộc có dạng bất đẳng thức thay vì là đẳng thức. Vậy : 


3) 18u® 6/(VJ) =1 nh) 


Người ta thường đưa bài toán này về bài toán vận tải (đóng) theo một trong hai trường hợp 
sau đây : 


1.Trường hợp thứ nhất là có tính đến sự thiệt hại bằng tiền khi thiếu một đơn vị hàng hoá 
Ở nút cầu j là rj (j=1 n) 


Lúc này người ta đưa thêm vào một nút cung giả (m+1) với nguồn cung là 
Sm+1l = S1—1đj — }i—1 8i 
và cước phí tương ứng là 
c(m+1) j = rj=1 n) 
Khi đó ta nhận được một bài toán vận tải (đóng) 
tH pin 6i 
S351 #ụ = đz(j = 1 — n) 


ng, Z¡ — 8;(# = 1 — m) 
2 0= 1.ev0m211J<e L<sĐ) 


ID 


2.Trường hợp thứ hai là không tính đến sự thiệt hại do thiếu hàng ở nút cầu 


Lúc này ta cũng đưa về bài toán vận tải (đóng) như trên, nhưng vì không tính đến sự thiệt 
hại nên mục tiêu sẽ là 


° Tn T 
mỉn ồ ;;_¡ È;;—¡ Cịij 
Ghi chú : 
Với bài toán vận tải mở, nguồn chuyển không hết sang các nhu cầu, người ta có thể tính 
thêm cước phí lưu kho ở mỗi nguồn cho mỗi đơn vị hàng là c¡ (n+1) (¡=1 m) . Hoàn toàn 


tương tự như trên, khi đưa bài toán này về bài toán vận tải (đóng) bằng cách thêm vào nút 
cầu giả (n+1) thì hàm mục tiêu trở thành 


: 1 
HH banh 3ê Cịjij 


Như vậy ta chỉ cần xét bài toán vận tải (đóng) 


š 1 
HH 327 3i E5 
» #j — s;¡( =l— mì) 


33-10 = đị(ƒ = 1 — n) 
#ọụ DU 1<ymjj}=1^šøn) 


t 


c- Bài toán vận tải có đường cấm 
Đây là bài toán vận tải nhưng không phải mỗi nguồn đều có cung nối với mọi đích. nghĩa là 
có đường cấm. Cách đưa về bài toán vận tải là dùng phương pháp M-lớn, tức là phương 
pháp phạt như sau : 
Gọi E là tập các cung không cấm, tức là các cung (¡,j), ¡ S, j D và bài toán có thêm điều kiện 
xij=0 với (1j) E 
ta đưa bài toán có các yêu cầu 
` m+1 m 
mỉn ồ”j¡ ồ „¡-¡ Cji 
và, Zj — s;( =l— mì) 


33-1 = đj(j = 1 —> n) 
#u 2 00=1->m.j=1->n) 
% — 0khi(ï,)) ⁄2 k (*) 


t 


về bài toán vận tải bằng cách đặt cước vận chuyển mới như sau : 


c¡nÕu(1,j) € E 
MnÕu({i,j) £ E 


đụ — { 


Ớ đây M là một số rất lớn, được coi là số lớn hơn mọi số gặp phải khi tính toán. 
Xét bài toán với cước phí mới như trên nhƯ sau : 
: m1 TH -= 
mỉn ồ”j¡ Ð „¡-¡ Ðji 
» Z¡ — s;( =l— mì) 


Le <=d,(0<1<9:B) 
œụ > 0(2=lL->mj]=1l>n) 


(3) 


t 


thì ta có : 
Định lý : 
Giả sử ø = [z¡]„.„là phương án vận chuyển tối ưu của (**) thì khi đó : 


1. Nếu z = 0V(i,j) £ E thì zlà phương án vận chuyển tối ưu của bài toán vận tải có 
đường cấm (*) 


2. Nếu tổn tại z„ị £ mà zạị > 0 thì bài toán vận tải có đường cấm (**) không có nhiệm 
chấp nhận được. 


d- Bài toán vận tải kèm chế biến trung gian 

Giả sử rằng trong mô hình vận tải có một số điểm nguồn, tức là điểm sản xuất, cho ra một 
số sản phẩm cần phải chế biến trước khi đến điểm cầu. Giả sử có =1 k điểm chế biến 
với khả năng chế biến là a đơn vị sản phẩm tương ứng. Gọi cước phí vận chuyển một 
đơn vị bán sản phẩm từ ¡ đến là c;; và chuyển một đơn vị sản phẩm từ đến j là c;:. Bài 


toán đặt ra là lập kế hoạch vận chuyển tất cả các sản phẩm qua chế biến đến tất cả các 
điểm cầu sao cho cước phí nhỏ nhất. 


Gọi xi j là lượng sản phẩm từ ¡ qua rồi qua j, ta cần tìm x=[ xi j]mkn sao cho : 


' m k n ? ?2 
min ồ,~¡ 3 ;A-i 2xi1CA ĐẸ CAj)®¡Aj 
k m › 
2 5A—1 3 2j—1 ý — Si( = 1 — mm) 
Tn k b 
3221 2 2= #àj = đị( = 1 — n) 
( 


3n bế" #j = đÀ(ÀA = 1 — È) 


t 


BÀI TOÁN DÒNG TRÊN MẠNG 


Mở đầu 

Nhiều bài toán quy hoạch tuyến tính có thể quy về bài toán làm cực tiểu phí tổn vận 
chuyển hàng trong một mạng (gồm các nút và các cung đường) sao cho đảm bảo được các 
nhu cầu ở một sỐ nút sau khi biết nguồn cung cấp tại một số nút khác. Các bài toán như 
vậy được gọi là các bài toán dòng trên mạng hay bài toán chuyển vận (TransShipment 
Problem). Đây là lớp bài toán quan trọng nhất và hay gặp nhất trong quy hoạch tuyển tính. 
Lớp này bao gồm các bài toán quen thuộc trong thực tế như : 

- Bài toán vận tải 

- Bài toán mạng điện 

- Bài toán mạng giao thông 

- Bài toán quản lý 

- Bài toán phân bổ vật tư 

- Bài toán bổ nhiệm 

- Bài toán kế hoạch tài chính 

- Bài toán đường ngắn nhất 


- Bài toán dòng lớn nhất 


Vì là một bài toán quy hoạch tuyến tính nên các bài toán dòng trên mạng có thể giải được 
bằng bất kỳ thuật toán nào giải được bài toán quy hoạch tuyến tính, chẳng hạn bằng thuật 
toán đơn hình như đã biết . Tuy nhiên, nếu tận dụng những cấu trúc đặc biệt của các bài 
toán dòng trên mạng sẽ làm cho phương pháp đơn hình đơn giản hơn và được thực hiện 
nhanh hơn. 

Phát biểu bài toán dòng trên mạng 


Mạng là một đồ thị có hướng ký hiệu G=(N,A), N là tập các nút, A là tập các cung, cùng 
một số thông tin về số lượng bổ sung như sau : 


. bi ( N) biểu thị nguồn từ ngoài vào nút ¡, gọi tắt là nguồn 

. uij biểu thị tải năng củỦa cung (¡,j) A 

. c¡j biểu thị cước phí cho một đơn vị của dòng trên cung (i,j) A 

. xij biểu thị lượng vận chuyển của dòng trên cung (¡,j) A 

Giá trị tuyệt đối |bi| được gọi là nhu cầu của nút ¡. Nếu bi>0 thì nút ¡ được gọi là điểm 
nguồn, nếu bi<0 thì nút ¡ được gọi là điểm hút. Một cách hoàn toàn tự nhiên người ta đặt hai 
điều kiện sau đây : 


a- Tổng lượng trên dòng vào nút ¡ bất kỳ phải bằng tổng lượng trên dòng ra khỏi nút ¡ (luật 
bảo toàn dòng). Như vậy : 


bị À3) ;crg) #äš = À3 2jeo() #u(¿ € N)Œ) 

Trong đó : 

I()= {nút j / cung (J,i) A} : những nút có cung nối đến nút i 

O()= {nút j / cung (i,j) A} : những nút có cung nối từ nút ¡ đến nó 

b- Dòng trên cung là không âm và không vượt quá tải năng của cung. Như vậy : 
0 <Szụ <S u¿VGjJ) € A(2) 


Mọi vectơ x có các thành phần xij, (¡,j) A, được gọi là một dòng. Dòng x thoả điều kiện (1) 
và (2) được gọi là dòng chấp nhận được. Lấy tổng của (1) theo các nút ¡ ta được : 


À ;cN b¿ = 0(3) 


Điều này có nghĩa là tổng dòng từ bên ngoài vào mạng phải bằng tổng dòng tỪ mạng ra 
ngoài. Nếu điều này điều này không thoả thì bài toán là không chấp nhận được. 


Mục tiêu của bài toán là làm cực tiểu cước phí dòng trên mạng, tức là : 


min „y4 cụ#i¡(4) 


trong đó cực tiểu lấy trên mọi dòng chấp nhận được. Như vậy ta nhận được một bài toán 
quy hoạch tuyến tính nhƯ sau : 
min Ð }\¡ ¡c4 Cjij 
bị Ð À3) jcrg) #ñ = 3 jco() #0? € N) 
0< Zj < œ¡¡V(1,j) cA 


QUY HOẠCH NGUYÊN 


Mở đầu 


Quy hoạch nguyên (Integer Programming), viết tắt là IP, là bài toán quy hoạch mà trong đó 
tất cả hoặc một phần các biến bị ràng buộc chỉ lấy giá trị nguyên. Trường hợp thứ nhất 
được gọi là quy hoạch nguyên hoàn toàn (Pure Integer Programming — PIP), trường hợp thứ 
hai được gọi là quy hoạch nguyên bộ phận (Mixed Integer Programming - MIP). Tuy vậy 
thuật ngữ ”quy hoạch nguyên” được dùng chung cho cả hai trường hợp. 

Mảng các bài toán có vẻ đơn giản nhất mà cũng là quan trọng nhất trong lớp các bài toán 
quy hoạch nguyên là các bài toán chọn các quyết định (chọn/không chọn). Chẳng hạn như 
bài toán bổ nhiệm, biến quyết định việc bổ nhiệm nhận giá trị như sau : 


1nÕu ngêi ¡ nhẼn c«ng viÖc j 
0nÕu ngêi ¡ kh«ng nhÈẼn c«ng viÖc j 
®ị — { 
Vì các biến quyết định thường chỉ nhận một trong hai giá trị nên bài toán này còn được gọi 
là bài toán quy hoạch nguyên nhị phân (Binary Integer Programming) . 


Một ý tưởng tự nhiên để giải bài toán quy hoạch nguyên là cứ giải như một bài toán quy 
hoạch tuyến tính tổng quát tạm bỏ qua ràng buộc biến phải nguyên. Khi tìm được phương 


án tối ưu thì sẽ làm tròn nó để được phương án tối ưu nguyên gần đúng. Phương pháp này 
có thể áp dụng trong thực tế nhưng phải chú ý đến hai nguy cơ sau đây : 


- Một là phương án tối ưu đã được làm tròn không chấp nhận được đối với bài toán quy 
hoạch nguyên. 


- Hai là phương án tối ưu đã được làm tròn chấp nhận được nhưng có thể giá trị mục tiêu 
tương ứng là rất xa với mục tiêu tối ưu của bài toán quy hoạch tuyến tính nguyên. 


Bài toán quy hoạch nguyên trong thực tế 
a- Bài toán balô 


Một nhà thám hiểm mang theo một balô chỉ chứa được một trọng lượng không quá b. Có n 
loại vật dụng phải mang theo. Mỗi vật loại vật ¡ có trọng lượng là ai và giá trị sử dụng là 
c¡. Hỏi ông ta phải chọn lựa các vật mang theo như thế nào để có giá trị sử dụng là lớn 
nhất ? 


Gọi xi (¡=1 n) là số lượng vật loại ¡ mà ông ta mang theo thì mô hình toán của bài toán balô 
này là quy hoạch nguyên nhƯ sau : 
2 
1n1aX2Z — Nai C¡#; 


334i 0ø <b 
z¡ > vh nguy°n(¿ =  — m) 


Về mặt toán học thì nếu hàm mục tiêu là min z hoặc ràng buộc là đẳng thức thì bài toán 
cũng gọi là bài toán balô. Bài toán balô có dạng đặc biệt và đơn giản vì chỉ có một ràng buộc 
ngoài ràng buộc dấu và tính nguyên. Người ta nghiên cứu được nhiều cách giải riêng cho bài 
toán và đưa bài toán quy hoạch nguyên về bài toán balô để giải. 


b- Bài toán sản xuất có lệ phí cố định 

Giả sử một nhà máy có kế hoạch sẽ sản xuất n sản phẩm. Chi phí sản xuất sản phẩm 

j=1 n gồm lệ phí cố định kj , không phụ thuộc vào số lượng sản phẩm j, và cước phí cj đối 
với mỗi đơn vị sản phẩm j. 


Gọi xj 0 là lượng sản phẩm j=1 n sẽ sản xuất thì chi phí sản xuất sản phẩm j sẽ là : 


k;+c;zxmDux; >0 
0nÕu Xj= 0 


c;(;) = { 


mục tiêu sản xuất với chi phí cực tiểu sẽ là : 
h .. T 
mỉnz = ồ},-¡ €/(#) 


Trong trường hợp này hàm mục tiêu z là hàm phi tuyến với các đối số là xj (=1 n) mặc dù 
các ràng buộc thực tế như nguyên liệu, thị truỜng,.... đều là tuyến tính nên bài toán rất khó 
giải. Người ta có thể đưa bài toán này về bài toán quy hoạch tuyến tính nguyên bộ phận 
bằng cách đưa vào các biến phụ nhị phân nhƯ sau : 


l1nÕu x;>0 
0nÕu xi= 0 
(1) 


‹® 

S, 
| 

¬ 


Để biểu thị yj (J=1 n) là biến nhị phân độc lập, không phụ thuộc vào xj như trong (1) người 
ta đưa vào một ràng buộc tuyến tính như sau : 


xj Myj =1 n) 


ở đây M>0 và rất lớn để ràng buỘc xj là thừa. Khi đó hàm mục tiêu và ràng buộc trên trở 
thành : 


mỉnz = 3};_¡(7 + c;Z;) 


0<z; < My, 
[0 


(2) 


Thật vậy : 


- Nếu xịj > 0 thì yj không thể bằng 0 nên yj =1 


- Nếu xịj = 0 thì yj = 0 hoặc yj=1 

Nhưng vì kj>0 ( nếu kj= 0 thì không cần đưa vào biến phụ yj) và hàm mục tiêu là min z nên 
Ở thuật toán tìm phương án tối ưu luôn lấy yj=0 vì phương án với xj=0 và yj=1 không thể là 
tối ưu. Khi viết đủ các ràng buộc tuyến tính khác vào ta được bài toán quy hoạch tuyến tính 
nguyên bộ phận. 


CÂU HỒI CHƯƠNG 4 

1- Trình bày chiến lược bị trội hơn. 

2- Trình bày chiến lược MaxiMin và MiniMax. 

3- Xây dựng quy hoạch tuyến tính trong trường hợp không có nghiệm ổn định. 
4- Trình bày các giai đoạn giải bài toán vận tải. 

BÀI TẬP CHƯƠNG 4 


1- Tìm phương án tối ưu cho bài toán lý thuyết trò chơi có ma trận điểm được cho như sau : 


2 3 -2 -1 
-1 Đ 4 -2 
-2 -ÐD 0 3 


2- Giải bài toán vận tải có ma trận cước phí 


60 70 40 30 
100 . 1 4 3 
30 Đ 3 2 6 


20 6 2 1 5 


Bài tập tổng hợp 


BÀI TẬP TỔNG HỢP 
I- Xí nghiệp sản xuất giấy có 3 phân xưởng. Do trang bị kỹ thuật khác nhau nên mức 


hao phí tre gỖ, axit để sản xuất một tấn giấy thành phẩm cũng khác nhau. Mức hao 
phí được cho trong bảng dưới đây : 


Mức hao phí nguyên liệu cho 1 tấn giấy 


Nguyên liệu P.Xưởng I P.Xưởng II P.Xưởng III 
Tre gỗ 1,4 (tấn) 1,3 1,7 
Axit 0,1 0,12 0,15 


Số lƯỢng tre gỖ có trong năm là 1.500.000 tấn, Axit là 100.000 tấn. 
Yêu cầu 


1. Xây dựng mô hình sao cho tổng sỐ giấy sản xuất trong năm của xí nghiệp là nhiều 
nhất. 


2. Xây dựng mô hình bài toán đối ngẫu với mô hình toán của câu 1. 


3. Tìm phương án tối ưu ứng với mô hình toán ở câu 1. TỪ đó suy ra số tấn giấy của 
mỗi phân xưởng cần sản xuất trong năm. 


4. Áp dụng kết quả bài toán đối ngẫu để từ bảng đơn hình tối Ưu câu 3 suy ra 
phương án tối ưu cho bài toán đối ngẫu câu 2. 


II- Một xí nghiệp có thể sản xuất bốn loại mặt hàng xuất khẩu H1, H2, H3, H4. Để 
sản xuất 4 loại mặt hàng này, xí nghiệp sử dụng 2 loại nguyên liệu N1, N2. Số 
nguyên liệu tối đa mà xí nghiệp huy động được tương ứng là 600kg và 800kg. Mức 
tiêu hao mỗi loại nguyên liệu để sản xuất một mặt hàng và lợi nhuận thu được được 
cho trong bảng sau : 


Định mức tiêu hao nguyên liệu và lợi nhuận H1 H2 H3 

NÑI1 0,5 0,2 0,3 

N2 0,1 0,4 0,2 

Lợi nhuận 0,8 0,3 0,5 
Yêu cầu 


1- Lập mô hình để xí nghiệp sản xuất đạt lợi nhuận cao nhất. 


2- Xây dựng bài toán đối ngẫu ứng với mô hình toán ở câu 1. 


H4 
0,4 
0,5 


0,4 


3- Áp dụng thuật toán đơn hình cải tiến và kết quả đối ngẫu để tìm các phương án 


tối ưu cho cả 2 mô hình. 


II- Xí nghiệp cơ khí Hùng VưƠng có 32 công nhân nam và 20 công nhân nỮ. Xí 


nghiệp có 2 loại máy : cắt và tiện. Năng suất trung bình của các công nhân đối với 


mỗi loại máy được cho trong bảng bên dưới đây : 


Năng suất công việc công nhần nam công nhân nữ 
Máy cắt 30 chỉ tiẾt/giỜ 22 chỉ tiẾt/giỜ 
Máy tiện 25 chỉ tiết/giỜ 20 chỉ tiết/giỜ 


Biết rằng trong ngày cắt được bao nhiêu chỉ tiết thì tiện hết bấy nhiêu chỉ tiết 


Yêu cầu 
1- Lập mô hình để xí nghiệp sản xuất được nhiều sản phẩm nhất. 


2- Lập mô hình đối ngẫu ứng với mô hình câu 1. 


3- Áp dụng thuật toán đơn hình cải tiến và kết quả đối ngẫu để tìm phương án tối 
ưu cho cả 2 mô hình toán trên. 


IV- Một công ty chuyên sản xuất 3 loại sản phẩm A, B, C. Trong đó nguyên liệu để 
sản xuất ra 3 loại sản phẩm trên được nhập về từ 2 nguồn N1, N2. Chi phí cho mỗi 
đơn vị nguyên liệu nhập từ nguồn N1 là 100000 USD và nguồn N2 là 90000 USD. 


Các loại sản phẩm sản xuất cần các đơn vị nguyên liệu của từng nguồn được cho 
trong bảng sau : 


Nguồn nguyên liệu Loại sản phẩm 

A B C 
NI 1000 2000 3000 
N2 2000 1000 2000 


Số lượng tối thiểu sản phẩm loại A cần sản xuất trong thời gian tới là 20000 , sản 
phẩm loại B là 18000, sản phẩm loại C là 15000. 


Yêu cầu 


1- Lập mô hình để tổng chi phí sản xuất mà công ty bỏ ra là nhỏ nhất mà vẫn đảm 
bảo yêu cầu về sản phẩm. 


2- Lập mô hình để công ty sản xuất đạt doanh thu cao nhất 
3. Tìm phương án tối ưu cho cả 2 mô hình. 


V- Một cơ sở dự định sản xuất tối đa trong một ngày 500 ổ bánh mì dài và 500 ổ bánh 
mì tròn, muốn đạt lợi nhuận nhiều nhất, với những điều kiện như sau : 


- Giá bán một ổ bánh mì dài làm từ 400 gam bột là 325 đồng, một ổ bánh mì tròn làm 
từ 250 gam bột là 220 đồng. 


- Số lượng bột được cung cấp tối đa trong ngày là 225 kg với giá mỗi kg là 300 đồng. 


- Lò nướng bánh cho phép nướng 75 ổ bánh mì dài hay 100 ổ bánh mì tròn trong mỘt 
giờ nhưng không thể nướng hai loại cùng một lúc. Lò nướng hoạt động tối đa 8 giờ 
trong một ngày. 


Yêu cầu 

1- Lập mô hình cho bài toán nêu trên. 

2- Xây dựng bài toán đối ngẫu cho bài toán trên. 

3- Tìm phương án tối ưu cho cả hai bài toán. 

VI- Ba xí nghiệp A, B, C cùng có thể sản xuất áo vét và quần. Khả năng sản xuất của 
mỗi xí nghiệp nhƯ sau : Khi đầu tư 1000USD vào xí nghiệp A thì thu được 35 áo vét 
và 45 quần ; vào xí nghiệp B thì thu được 40 áo vét và 42 quần ; vào xí nghiệp C thì 


thu được 43 áo vét và 30 quần. LƯợng vải và và giờ công để sản xuất được cho trong 
bảng sau : 


Xí 
nghiệ 
phiệp A B C 
Vải/Giờ Vải/Giờ Vải/Giờ 
1áo s& 16 18 
vết 3,5m 20 giỜ 4m giờ 3,8m giờ 
A vx 12 15 
1 quần 2,8m 10 gIỜ 2,6m giờ 2,om giỜ 


Tổng số vải huy động được là 10000m. 


Tổng số giờ công huy động được là 52000 giỜ. 

Theo hợp đồng thì tối thiểu phải có 1500 bộ quần áo, nếu lẻ bộ thì quần là dễ bán. 
Hãy lập kế hoạch đầu tư vào mỗi xí nghiệp bao nhiêu vốn để : 

- Hoàn thành hợp đồng 

- Không khó khăn về tiêu thụ 

- Không bị động về vải và giờ công 


- Tổng số vốn đầu tư là nhỏ nhất 


